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Abstract 

In paper a generalized definition of superpotential has proposed, which connects two 
one-dimensional potentials Vi and V2 with discrete energy spectra completely and where: 

• for definition of an energy of factorization an arbitrary level of the energy spectrum 
of the given Vi is used and a function of factorization is defined concerning a bound 
(ground or excited) state at this energy level, 

• for the definition of the energy of factorization an arbitrary energy (which can be not 
coincident with levels of the spectrum of Vi) is used and the function of factorization 
is defined concerning an unbound (or non-normalizable) state at this energy. 

It has shown, that for the unknown superpotential such its definition follows directly from 
a solution of Riccati equation at the given potential Vi . 

Using the arbitrary bound state in the construction of the superpotential, the SUSY 
QM methods at their effectiveness in construction of new exactly solvable potentials on 
the basis of the given potential Vi , in detail calculations of spectral characteristics under 
the deformation of Vi have been coming to a level of methods of inverse problem. So, if 
as the given potential Vi to choose a rectangular well with finite width and infinitely high 
walls, then we reconstruct by SUSY QM approach all pictures of deformation (without 
displacement of the levels in the spectrum) of this potential and its wave functions of the 
lowest bound states, which were obtained early by the methods of the inverse problem. 
Here, interdependence between parameters of the deformation for the methods of SUSY 
QM and the inverse problem has found, an analysis of a behavior of wave functions and the 
potential under the deformation has fulfilled, a classification has proposed for zero-points 
of the potential, nodes of the deformed wave functions, points, where wave functions are 
not deformed, an analysis of angles of wave functions leaving from such points has fulfilled. 

Using the unbound and non-normalizable states at the arbitrary energy of factorization 
in the definition of the superpotential, we obtain new types of deformations. So, using 
only one superpotential, one can join two potentials, which have the real energy spectra 
with own bound states and without coincident levels. 
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Designations 



In paper we shall use the following designations. 

• The first potential T4, on the basis of which we construct new potentials using SUSY- 
transformations, we shall name as the starting potential with number "1 ". 

• (pf^^ is wave function (WF or WFs — in the plural case, partial or general solution) of 
the potential T4 with the number k, where instead of index i at the bottom we shall use 
one index from the following: 

— n — as the number of the level of discrete energy spectrum of the potential with number 
k, if this level does not coincide with energy of factorization £ (i. e. the energy, at which 
the superpotential is defined); 

— m — as the number of the level of the discrete energy spectrum of the potential Vfe with 
the number k, if this level coincides with the energy of factorization £; 

— w — as the index at energy, which coincides with the energy of factorization £ and does 
not coincide with any level of the energy spectrum of the starting potential Vi. 

• Ef''' is the level with the number i (corresponding to WF f^^^^) of the energy spectrum 
for the potential Vk- We shall number the levels of the discrete spectrum so: 

— the lowest level, concerned with the ground bound state, by number "1"; 

— the next levels, located higher and described excited bound states, by numbers "2" and 



• Wl ' is the superpotential, for which as the energy of factorization we shall choose the 
energy with index i, and as a function of factorization wc shall choose WF ^"i' (of bound 
or unbound state) of the potential Vk- If the function W^/^'' will be considered with only 
one number k in the text, then we omit this number sometimes. 

• We shall choose such definition of operators Ai and Af: 



• ipi^^'> is WF of unbound or non-normalizable state, which we shall denote by stroke on 
top. 



so on. 




(1) 
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1 Introduction 



Supersymmetric theories are ones from the most modern and vigorously developed topics of 
physics of particles and interactions. Usually, under a word "super symmetry" we have in 
mind a relation established between bosons and fermions, which allows to consider the boson 
and the fermion as two different manifestations of an unified particle [1]. It is interesting 
to note, that according to one of main principles of quantum mechanics (without applying 
a formalism of field theories), a conception of "wave" (which we relate to interactions in the 
simplest understanding) and a conception of "particle" (under which we keep in mind the matter 
in the simplest understanding) can be considered as two different manifestations of only one 
object — "particle- wave" . Mathematically, for description of the particles and interactions the 
supersymmetric models (for example, such as Minimal Supersymmetric Standard Model [2, 3]), 
based on a modern apparatus of quantum field theory mainly with application of algebras of 
supersymmetry [4], are used. It makes study of properties of interdependence between fields 
of the matter (as the fermion fields) and fields of the interactions (as the boson fields) enough 
complicated. 

In 1981 E. Witten in his famous paper [5] constructed a simple example of one-dimensional 
supersymmetric system, which has no any reference to the theory of fields at all! It had opened 
a possibility to study absolutely independently properties of supersymmetric transformations 
between partners on the basis of the simplest examples, without using formalism of the field 
theories (and without formalism for description of spin) . Since then, after publishing of some 
fundamental papers, a new topic of physics had been created and has been developed vigorously, 
which is named as supersymmetric quantum mechanics (SUSY QM) [6]. 

An essential interest in SUSY QM can be explained by that it allows to study in details 
the properties of the supersymmetric transformations, using the simplest quantum systems. Of 
cause, then such results can be applied in the field theories up to cosmology [7] and quantum 
gravity. So, in the paper [8] (see p. 2917-2918) C. V. Sukumar clearly explained such connection 
so (also see [9, 10] of other authors): "It is well known that d = 1 quantum mechanics is 
formally equivalent to the d = 1 quantum field theory with the identification x ^ if, p ^ dt^p 
and canonical quantisation of the field ip leads to the usual commutation relation between x and 
p..." 

Prom another side, the second important perspective of SUSY QM has been opened: it 
improves essentially methods of quantum mechanics directly. Keeping in the frameworks of 

quantum mechanics, SUSY QM proposes to researcher its own simple mathematically and 
clear enough methods for calculation of spectral characteristics of quantum systems, it gives 
answers to questions and vagueness, putted early by methods of direct and inverse problems of 
quantum mechanics. First of all, one note a method of factorization (principally, introduced at 
the first time by A. A. Andrianov, N. V. Borisov and M. V. loffe in [11] and independently by 
C. V. Sukumar at first in [12] and then in [8, 13, 14, 15, 16, 17] where a main formalism of such a 
method was constructed, however, the algorithm of factorization, started else by P. A. M. Dirac 
in [18] and E. Schrodinger in [19], can be found earlier for solution of the spectral problem with 
some simplest one-dimensional potentials in quantum mechanics, for example, see [20, 21]), 
formulation of which has been causing a large number of papers in SUSY QM, search of new 
methods and opening of new types of exactly solvable potentials on the basis of intensive study of 
Darboux transformations (first of all, one can point out papers of B. F. Samsonov, V. G. Bagrov 
with their colleagues at development of a general formalism of these transformations with their 
application in a number of different problems [22, 23, 24, 25, 26], it is necessary to note papers 
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of other researchers at their investigations of these transformations with apphcation to theory 
of solitons [27, 28, 29], in approaches of nonlinear supersymmetry [30, 31, 32], in solution 
of matrix Schrddinger equations [33], in solution of system of coupled discrete Schrddinger 
equations [34] , in generalization of the formalism of these transformations for three-dimensional 
and many- dimensional spaces [11, 35, 36, 37], in scattering theory [38, 23, 25], in improvement 
of the methods of the inverse problem on the basis of these transformations [39] , in study of 
interdependence between these transformations, the method of factorization and the methods 
of the inverse problem [35, 40, 41], in intriguing and having large perspectives non- stationary 
approaches [42, 43], one can note approach of the inverse Darhoux transformations in description 
of some types of exactly solvable potentials [44] , and also interesting review about such Darboux 
transformations [45] and historical paper [46]). 

An essential progress was achieved early in development of the methods of the inverse 
problem in construction of the new exactly solvable potentials on the basis of known ones (see 
[47, 48, 49]). A vigorous accumulation of information about new types of the potentials with 
spectral characteristics, which have explicit analytical form, little by little forms an unified 
theory of exactly solvable models (or potentials) (see [47], p. 915). The methods of SUSY QM, 
after their opening of original solutions for new types of such potentials (for example, shape 
invariant potentials or SIP-potentials [50, 6] with translations of parameters [51, 52, 53], with 
scaling of the parameters [54, 55], with rotations of parameters [56], class of the potentials of 
Shahat and Spiridonov [57, 58, 51], which includes soliton-like reflectionless solutions without 
tunneling, search of the reflectionless potentials with barriers [59, 60, 61] (where an effect of 
reflectionless tunneling can appear, existence of which is discussed up to now), the potentials 
in approach of nonlinear supersymmetry [31, 32]), have proved their right to originality, and 
their effectiveness in detail calculations of spectral characteristics of the potentials comes little 
by little to a level of the methods of the inverse problem. During last two decades, one can see 
a continuous competition between the methods of the inverse problem and SUSY QM methods, 
which supplement each other [13, 15, 17]. However, perhaps, increase of number of the new 
SUSY QM methods looks more intensive, that one can explain by their comparative simplicity 
and availability. 

These arguments emphasize the interest in the further development of the new SUSY QM 
methods, which give to researcher the enlargcr possibilities in construction of the new exactly 
solvable potentials (on the basis of known potentials) in comparison with early known methods. 

In this paper we shall restrict ourselves by consideration of one-dimensional potentials with 
discrete spectra of energy completely. Here, we shall analyze more generalized definition of a 
superpotential, which we construct in the region of real values only and which connects two 
potentials SUSY-partncrs in frameworks of the most widely known approach of SUSY QM (for 
example, as in formalism [8, 6]), adding an arbitrariness into a choice of a bound state at an 
energy of factorization, which coincides with the ground or excited level of the energy spectrum 
of the given potential Vi, and introducing an arbitrary unbound and non-normalizable states at 
the arbitrary energy of factorization, which can be not coincident with the levels of the energy 
spectrum of Vi. We note, that similar generalizations on the construction of the superpotential 
were considered early (see [15] — more for the potentials with continuous energy spectra; [62] 
p. 378-386 — for the superpotential and the energy of factorization in the region of complex 
values). However, we see, that this question for the potentials with discrete energy spectra has 
not been studied yet practically, and our analysis has shown an extreme availability of taking 
into account of different types of states at the arbitrary energy of factorization in the definition 
of the superpotential, even if to find it in the region of real values only. 
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It turns out, that such generahzation follows directly from a form of the superpotential, if to 
find it as an exact solution of Riccati equation at a given potential Vi, own independent variant 
of which we present in Sec. 2 (see also original papers [8, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72] and 
monographs [73, 74]). It is necessary to note, that the Riccati equation admits an arbitrariness 
in a choice of boundary conditions in the definition of the superpotential, while the standard 
approach to the construction of the superpotential uses wave function of the bound (mainly, 
ground) state only and, therefore, it is partial. 

The generalization for definition of the superpotential requires introduction of the unbound 
or non-normalizahle states at the arbitrary energy of factorization, and a formalism of such 
states we present in Sec. 4 (also sec [8] p. 2925-2934, [6] p. 324-325, [49] p. 286-290, 293- 
297, 309-318, [62] p. 378-382, 382-389). In Sec. 5 we improve a method of construction of 
n-parametric family of isospectral potentials, developed early in [6] (see p. 326-328). 

Generalizing the definition of the superpotential on the basis of wave function of the arbi- 
trary excited bound state of the given potential Vi, we have achieved practically full analogous 
between SUSY QM methods and the methods of the inverse problem. So, in Sec. 6 we recon- 
struct all pictures of the deformation (without of displacement of the levels) of a rectangular 
well with finite width and infinitely high walls and its wave functions of the lowest bound states, 
which were obtained early in the review [47] (see p. 916-919) on the basis of the methods of the 
inverse problem. Here, we find interdependence between main parameters of deformation for 
the SUSY QM methods and the methods of the inverse problem, find answers (confirmed by 
analytical expressions) on some questions, putted early in this review. These arguments make 
the SUSY QM methods practically such effective as the methods of the inverse problem. 

We see, that use of the unbound states in the definition of the superpotential (in the energy 
region of the real values) opens new original perspectives for the SUSY QM methods. It 
turns out, that if to do not restrict oneself by the bound states only in the definition of the 
superpotential, but to construct it on the basis of a function, which is a general (or partial) 
solution (without implying the boundary conditions for bound states on it) of the Schrodinger 
equation with the starting potential Vi at the arbitrary energy (which can be not coincident 
with the energy levels of Vi), then we obtain new exactly solvable potentials which have own 
energy spectra with the bound states. In particular, using only one superpotential, one can 
join two potentials, energy spectra of which (with own bound states) have no one coincident 
level (it has found at the first time in the supersymmetric approach)! As a demonstration, 
in Sec. 7 we present two new approaches of construction of the new potentials, where as the 
starting potential Vi we use the rectangular well with infinitely high external walls. This opens 
a possibility to deform at once whole energy spectrum for the given potential (perhaps, by 
a given rule, for example, as for SIP-potcntials with scaling of the parameters). Perhaps, 
such a way will open a possibility to connect together boson and fermion components with 
absolutely different energy spectra in the supersymmetric quantum theory of fields, using only 
one superpotential. 

2 An exact analytical solution of a superpotential at the 
given potential 

To clarify, which the most general form an unknown superpotential can have, when we know 
only a form of the potential, concerning with this superpotential in the most prevailing for- 
mahsm (for example, as in [6] p. 275-277, [8] p. 2922-2923, 2925-2927, [11] p. 19-20, at an 
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energy of factorization, located not higher then the lowest level of energy spectrum of the given 
potential), we inevitably come to a. problem of a solution of Riccati equation, where as unknown 
function the superpotential is used. 



2.1 The superpotential as the solution of the Riccati equation 

Let's consider a quantum system with a potential Vi{x) and a discrete energy spectrum com- 
pletely. For such a system one can write down the Schrodinger equation: 



hMr^K^) = {-^,+Vi{x)yg\x)^E(^^^g\x), (2) 

where hi is Hamiltonian of the system, ipg'' (x) is wave function (WF) of a state with a number 
m of this system, E^^ is energy level, corresponding to the state with the number m, e[^^ is 
the lowest level of the energy spectrum with the number 1. For the given potential one can 
introduce a superpotential W{x), defining it by such a condition (for example, according to [6], 
p. 287-289): 

V,(.) = W'(.)-^ + E?K (3) 

Let the potential Vi{x) be exactly solvable, i. e. one can write WFs for all states and the 
energy spectrum of the system with such a potential in the explicit analytical form. Let's 
assume, that we know a form of the potential ^1(2;), its WFs and the energy spectrum. Also 
we assume, that wc do not know the superpotential W{x) and we shall find it. The condition 
(3), defining the unknown superpotential W{x), is the Riccati equation (see in the Introduction 
citations on some original papers and monographs). Let's find the form of the superpotential, 
solving this equation. 

We fulfill a substitution of variable x ^ u — u{x), defining it so: 

/W^^, (4) 
where on the function f{x) we impose arbitrariness in its choice. We obtain: 

du{x)^fdx, ^^f^. (5) 
dx du 

Then one can rewrite the equation (3) through new variables: 

V,{u)-E['> = WHu)-f(u)^. (6) 



Now let's consider a function of a from: 

a 

Wr{u) = { 



U — Un 

a 



U + Uq' 



at u <0, 

(7) 

at u>Q, 



where a > and Uq > are constant real parameters. By direct substitution one can make 
sure, that at a = 1 the potential 



VM = WUu)-'^ (8) 
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becomes zero (according to [59, 60], the function of the form (7) is only one from five possible 
solutions for the superpotential Wr, on the basis of which one can construct zero potential (8); 
therefore, the choice (7) further will give us only a partial solution in the definition of the su- 
perpotential W{x); therefore, for example, substitutions (A3) and y = l/<^ in [8], perhaps, lead 
to the partial definition of the superpotential U in the form (A12) at the energy of factorization 
E located not higher then the lowest level of the energy spectrum of the potential V) . In the 
definition of the function Wji{u) in the form (7) its divergence at point = is excluded by in- 
troduction of the non-zero parameter uq] however, the function Wii{u) has discontinuity at zero, 
nevertheless constructed on its basis the fimction Vr(-u) is continuous at this point, and, there- 
fore, on all axis u (for example, in contrast to the function y{x) in limit ^p{x) ^ at substitution 
y = l/(p in [8], p. 2934; exclusion of such a divergence for the function y{x) and/or (p{x) at 
passing of the variable x through the extreme value a; = can be interesting, because it allows 
to consider oscillations of the functions y{x) and (p{x); it is important to study a question 
of passing of these functions through their zero values or "nodes", nevertheless it turns out, 
that it is possible to define the superpotential at the energy of factorization, which coincides 
not only with the lowest level of the energy spectrum for the given potential or below it, but 
and higher then such level). Here, in future investigations, it is interesting to use a formal- 
ism for obtaining solutions for "the regulized one- dimensional form of the centrifugal potential 
V{x) = a/{x + ieY" where a is a real strength and e a real constant that removes the singular- 
ity at the origin, proposed in [75] (see p. 32-34), as a possible further complex extensions of a 
definition of the function Wr{u) in search of new types of exact solutions of Riccati equation 
by the approach proposed in this section. 

For simplicity of analysis further we shall consider the function Vr{u) on the positive semi- 
axis w > 0. Subtracting this function from the expression (6), we obtain: 

VM - Ef' = WHU) - IXuf^ - nu){wl(n) - «), (9) 

here multiplication of Vr{u) on square of the function f{u) leaves this function as zero. 
Now fix f{u) by such a condition: 

f{u)WR{u) = W{x). (10) 

Then from (9) we find: 

- E? = P(u)wi(u) - /(«)^ - f{u)wUu) + fHu)"!^ = 

du du 1 du 



Prom here we find f{u): 

df\u) _ 2{Ei'^ -V,{x)) 



:i2) 



du Wr{u) 

Thus, we have found the function f{u) in dependence on the known potential Vi{x) and the 
superpotential Wr{u). However, it is failed to integrate explicitly this expression else, because 
implicit dependence of the function Wr{u) on the variable / is included into the right part of 
(12), and Vi{x) depends on x. 
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Taking into account (5), we rewrite (12) by such a way: 

W^i^(«)^^ = 2Wn{u)f{u)^f{u) = '^Wniu)^ = 2{eP - V,{x 

Wn{u)^ = E['^-V,{x). (13) 

We have obtained the equation, in which the left and right parts depend only on one variable 
u or X. Further, using the solution (7) for non-zero superpotential Wr (at m > 0), we obtain: 

a df{u) a d du{x) a d^u{x) (i) 



u + Uq dx u + Uodx dx u + Uq dx"^ 

or (at a — 1) 

+ V,{x)j (u{x) + uo) = e['^ (u{x) + uo). (14) 

Thus, we have obtained the Schrodingcr equation for the potential Vi{x)\ 

Solving the Schrodinger equation at selected energy, one can find a general form of the 
function u{x) + uq. Usually, for potentials with infinitely high external walls the bound states 
cause a physical interest, WFs of which tend to zero in asymptotic space regions. Such boundary 
condition, imposed on WFs, introduces restrictions on possible energy values, defining the 
discrete spectrum. 

The level e[^^ is the lowest level of the energy spectrum for the potential Vi{x). Therefore, 
this level is eigenvalue of the Hamiltonian H^^\ for which there is non-zero solution for the 
function u{x) + Uq, being the WF. According to set of the problem, the value e[^^ and WF 
concerned with it are known to us. One can write: 

u(x) = Coipi\x)-Uo, 

^ dujx) _ ^ d^^^\x) 
dx ^ dx ^ _^ 



[X 



Co = const. 



Note, that the found solution for the superpotential Wn{x) has obtained only in such interval 
of values x, where u{x) > {u{x) is the continuous function in this interval), and therefore such 
a determination of the superpotential does not contain divergences. It agrees with a behavior 
of WF of the bound state at the lowest energy level, which has no nodes. 
One can define a form of a potential-partner V2{x) so: 

V,{x) = W\x) + '^^ + E^^\ (16) 

Usually, if the superpotential W{w) is defined in the form of fourth expression from (15), 
then the energy level E^"^ is named as an energy of factorization (let's denote it as £^), and 
the function i^'i\x), defining the superpotential by such a way, is named as a function of 
factorization. 
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Conclusion: 

The solution for the superpotential was obtained by resolving the Riccati equa- 
tion. In obtaining of the equation (I4) for determination of the superpotential we 

do not use the boundary conditions, which can be imposed on the function u{x). Further, ap- 
plying the boundary conditions, we define a state, described by the function u{x), as bound one 
at the energy of factorization, coincident with the lowest energy level e'^\ Therefore, the so- 
lution (15), constructed on the basis of WF of the bound state at the energy of factorization, 
coincident with the lowest level e'^\ can be a partial solution in definition of the superpotential 
W{x). 



2.2 Arbitrariness in a choice of the energy of factorization 

Now let's assume, that for considered above quantum system with the potential Vi{x) one can 
define the superpotential W{x) instead of the condition (3) by the following: 

V,{x) = W\x) - + ) + AE^'\ (17) 

where Ai?*^^^ is arbitrary real number. 

Assuming WF and energy spectrum for the potential Vi as known to us, we find the su- 
perpotential in this case. Then in the approach, presented in the previous section, we obtain 
instead of (14) the following equation: 



^ + T4(x)) (u{x) + u^) = ) + A£;(^)) (u{x) + t^o). (18) 

We again obtain the Schrodinger equation with the potential V\{x). From (18) we find: 

• If the energy Ei''^ + A£'*^^) coincides with arbitrary level E!^^ with a number m of the 
energy spectrum for the potential Vi, then it is eigenvalue of the Hamiltonian H^^\ for 
which there is non-zero solution u{x) + uq of the equation (18), described a bound state 
for Vi at the selected value AE^'^K Then the solution for the superpotential iy(a;) exists, 
where Ei'^ + AE^^^^ is the energy of factorization: 

_ du{x) _ ^ dipg^x) 



dx ^ ^ dx 

= -swT^=-^::^' (19) 



[X 



Wix) = f(:c)Wniu{x))^-^^^^^^ = -^log\^g\x)\, 



Cm — const. 



^^^{x) dx dx 



• If the value E^ + AE'^'^'> does not coincide with any level E':^ from the energy spectrum 
of the potential Vi, then it cannot be eigenvalue of the Hamiltonian hi and, therefore, 
there is no any non-zero solution u{x) + Uq of the equation (18), described the bound 
state for Vi at the selected value AE'^^^ (with interdependence (17) between the potential 
Vi{x) and the superpotential W{x)). 
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• Each separate known solution for wave function with its level of the energy spectrum 
of the potential Vi gives a new independent solution for the superpotential W{x), and 
concerned with it the energy of factorization. 

However, so defined superpotential can have divergent values at points-nodes of WFs of 
exited bound states for Vi, it turns out, that by use of it one can construct new exactly solvable 
potentials with bound states and without divergences. Therefore, further in this paper we shall 
suppose a possibility of definition of the superpotential with divergences at isolated points. 

2.3 An exact solution for the superpotential with next number in 
the hierarchy 

Let's consider three Hamiltonians and /in+2 for three potentials Vn{x), Vn+i{x) and 

Vn+2{x) with discrete energy spectra completely: 

k+ivt^'\x) = + = E^^+^)^t+^){xl (20) 

K^2^t'\x) = (-^ + K+2(x))^r^)(x) = Er''^t'\xl 



where (x), ip^^'^^^x) and ip"^^'^\x) are WFs for states with numbers /, m and k for three 
Hamiltonians; E'f-\ E^+^^ and E^"^^^ are levels of the energy spectra of these Hamiltonians. 
Let's assume, that the Hamiltonians hn and hn+i connected between themselves by a superpo- 
tential W^'^\x), and the Hamiltonians hn+i and hn+2 — by a superpotential W'^'^'^^^x), where 
the superpotentials are defined so: 

where and are the lowest levels of the energy spectra. Prom such definition we 

obtain: 

From (21) and (22) one can find equation for connection the superpotentials W^^'^ and 
with neighboring numbers n and n + 1 in the hierarchy: 

(„....(,,)^ + . _ ''^'";;(-) + (23) 

where kn+i — e["'^^^ — e[^\ This equation is widely used in different topics of SUSY QM, and, 
in particular, in parasupersymmetric quantum mechanics or PSUSY (one can see a brief review 
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of PSUSY with a literature list in [6] p. 370-377, at fist time PSUSY of order 2 was introduced 
in [76], PSUSY of arbitrary order p was proposed in [77, 78]). 

Let's we know the superpotential W'^'^^x) and WF of the potential Vn{x) at the level, 
coincident with the energy of factorization of the second superpotential W^'^\x). Then as 
it was shown in the previous section, the equation (23) can be resolved exactly analytically 
concerning the unknown superpotential W^"'~^^^ with the next number n + 1. We can find 
exactly analytically as many new solutions for the superpotential W^'^'^^\ as many WFs with 
levels corresponding to them ( and not coincident with the energy of factorization of the first 
superpotential W^'^'') for the potential Vn{x) we know (it improves the early known approach for 
determination of exact solutions for the superpotential with the next number in the hierarchy 
on the basis of known superpotential with the previous number). So, if the potential Vn{x) is 
exactly solvable and we know all its spectral characteristics, then the method in sec. 2 gives a 
whole set of the new exact solutions for the superpotential W^'^'^^\ Partially, such approach can 
be used for construction of new exactly solvable models in PSUSY QM. 

3 A variety of potentials-partners V2 to one potential Vi 



Arbitrariness in the choice of the energy of factorization gives a set of different potentials V2 
partners to one potential Vi. Let's consider this question in more details, when superpotential 
is constructed on the basis of wave function of an arbitrary bound state of Vi. For the given 
Vi we have a Hamiltonian hi of a form: 

hiV^:\x) = (-^, + Vi{x))^l^\x) = Ei\(:\x). (24) 



dx"^ 

According to analysis in the previous section, for the same Vi one can construct the superpo- 
tential W by different ways: 

Vi{x) = W\x) - ^ + E? = W\x) - ^ + + A£;«. (25) 



Including operators: 



A+^—^ + W(x), A+ ^—^ + W(x), 
dx dx 



(26) 



one can write: 

k^(^\x) = (a^a + e^%^P^ 



X) 



A+A + e['^ + AEW]ip(^\x) = E^'^ip'^^x 



(27) 



Using the found W and W, one can define new V2 and V2, which are partners to the same Vi. 

dx (-28) 

V2{x) = W\x) + + Ei'^ + AE^'\ 

dx 
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Conclusion: 

For one given potential Vi one can define a whole set of superpotentials W ( equals to number of 
known energy levels in the spectrum for Vi), each of which defines the new potential V2, which 
is the partner to the same Vi. 

3.1 Whether are the potentials V2 cind V2 exactly solvable and what 
are their spectral characteristics? 

Let's find out, whether the potentials V2 and V2, constructed as the partners to the same Vi by 
use of two different superpotentials W and W, are exactly solvable. 

For analysis of the potential V2, constructed at the energy of factorization S = Ei'\ we use 
expression for the Hamiltonian hi through operators A and A'^: 

h^^^\x) = (ma + E?y^\x) = (29) 

Let's act on this equation by the operator A from the left: 

A(hiif^^){x)) = a[a+a + ))<^W(x) = 

= {aA+^ (x)) + E? (Av^W (x)) = (30) 

= (aA-^ + E«) (M^)(a:)) = /.2(M^)(^)) = Eg\Aip^^){ 



Prom (30) one can see, that values E^^ represent eigenvalues of the Hamiltonian /i2, and 
functions Aip^^'^ represent eigenf unctions of this operator to normalizing constant. It proves, 
that new potential V2 is exactly solvable. We write down wave functions and the energy 
spectrum of this potential: 

V2{x)^W\x) + ^^ + E?, Ei^)^Ei^\ ^3^) 
(p^^\x) = CmA(p^^\x), Cm = const. 

For the analysis of the potential V2 at the energy of factorization S = e[^^ + AE^^\ we use 
expression for the Hamiltonian hi through operators A and A'^: 

hM:\x) = (A+A + e[''> + AE(''>y^\x) = E(^^ip(^\x). (32) 

Acting on this equation from the left by the operator A, one can obtain the Schrodinger equation 
again. Prom here we find, that new potential V2 is exactly solvable also, and its wave functions 
<^^^ and energy spectrum E^^ have forms: 

V2ix) = W\x) + ^ + Ei'^ + AE('\ Ei) ^ Ei^\ (33) 
(^(^) (x) = C^^c^W (x) , Cm = const. 

Conclusion: 

The set of potentials V2 is exactly solvable completely. Irrespective of the arbitrariness in the 
choice of the energy of factorization, all potentials V2 have coincided energy spectra (with possible 
exception of selected levels), and their wave functions are determined explicitly from (31) or 
(33). 
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3.2 Absence of additional levels in the spectra of the potentials V2 
and V2, and exception of a level coincided with the energy of 
factorization 

In the beginning we shall analyze wave functions for levels of the energy spectra of the potentials 
V2 and V2 in the case when these levels do not coincide with the energy of factorization S, i. e. 
at E^^ i^E = ^ for and at E^^^ + A£;(i) for V2. According to (31) and (33), 

for each level E^'^ of potential Vi there is as a minimum one non-zero solution for wave function 
at the level E^ , which coincides with for the potential V2 and there is as a minimum 
one non-zero solution for wave function at the level E^jP , which coincides with E!^\ for the 
potential V2. Acting by the operators and on the equations: 

h,^f{x) = + E?y^^\x) = E^^^fix), ^^^^ 

we find, that for each level E!^^ for V2 and for each level E^'' for V2 there is as a minimum one 
non-zero solution for wave function for Vi at the level which coincides with them. Thus, 
we come to exact one-to-one correspondence between all levels (with non-zero wave functions) 

of the spectra of energy for the potentials Vi, V2 and V2 (i. e. there are no any additional 
level for V2 and V2, which do not exist in the spectrum of Vi, at the chosen definition of the 
superpotential) : 

= Ei^^ = ^ 8. (35) 

Now let's consider the levels of the potentials V2 and V2, which coincide with the energy of 
factorization S, i. e. the level Ej^,} = 8 for V2 aX 8 = Ei''' and the level Ef^ = 8 for V2 at 
8 — e[^^ -\- AE^^\ Analyzing all possible states at these levels, we conclude: 

• The Hamiltonian h2 cannot have a state with non-zero wave function at the energy level 
EI^\ which coincides with the energy of factorization 8 = e[^\ Therefore, the energy 
spectrum of the potential V2 coincides completely with the energy spectrum of the potential 
V\ with possible exception of the single level E^'' . 

• The Hamiltonian /12 cannot have a state with non-zero wave function at the energy level 
E^^ , which coincides with the energy of factorization 8 — E^^ -\- AE'^^^ . Therefore, the 
energy spectrum of the potential V2 coincides completely with the energy spectrum of the 
potential V\ with possible exception of the singe level E[^^ + AE(^\ 

3.3 How do shapes of the potentials V2 and V2 differ? 

Let's find out, how do the shapes of the potentials V2 and V2 differ. Taking into account (25) 
and (28), we obtain: 

y, = V,-2^ = V-,-2^f. (36) 

Prom here find: 

V2 = V2 + 2^{W -W). (37) 
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According to (15) and (19), we see, that the superpotentials W and W are expressed through 
different wave functions and, therefore, the diflFerence W — W is not reduced to constant. One 
can write down: 

V2 = V2 + 2— log 



'■Pm 



(3J 



Conclusion: 

The potential V2 at its shape does not coincide with the potential V2, and, therefore, it is the 
new exactly solvable potential. 



4 New unbound and non-normalizable states in discrete 
spectra 

4.1 A partial solution of wave function of the unbound state 

According to the approach from sec. 2.1, the superpotential W at the given potential Vi can 
be constructed at arbitrary energy of factorization 8. If the energy of factorization coincides 
with arbitrary selected level Ejf^ of the bound state of Vi (as in sec. 3), then we write down: 

W^- = W^i'^ = -£log|(^«|, (39) 

where </?^^ is WF of the bound state with number m for Vi and the index m is added to the 
superpotential. The superpotential is connected with the potential V\ so: 

and a new potential V2 can be defined so: 

According to sec. 3.2, the potential V2 has no bound states at the level E^j^K However, such 
a function: 

^-J'^'^ - 4) (42) 

is a solution of the Schrodinger equation for this potential at this level. The wave function 
describes the bound state for Vi, tending to zero in both directions in asymptotic regions. 
Therefore, the function i^rrS^^ tends to its infinite values in asymptotic limits and can not 
describe a bound state for V2. Usually, such states are not used in practice, because they do 
not contain obvious physical sense for potentials with discrete energy spectra. However, it turns 
out, that on the basis of such states one can improve essentially earlier known supersymmetric 
methods of construction of new exactly solvable potentials at the given potential Vi, building 
new potentials with own bound states. Therefore, further we shall name such states as unbound 
(or non-normalizable) states in the discrete energy spectrum, and the functions, describing them, 
as wave functions of such states. If the function ^p^^ has nodes, then the function ^PrrS^'' has 
divergences at coordinates of these nodes. In such a case one can name the state, described by 
the function (pm^'^\ as discontinuous (broken) state. 
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From (39) and (42) one can write down: 



= llog\^^^'% (43) 



Conclusion: 

SUSY-algorithm of construction of new potential V2 with the neoct number in hierarchy with 
the given potential Vi destroys one bound state at the level, which coincides with the energy of 
factorization, and introduces a new unbound or discontinuous state (both partial, and general 
solutions for its WF) at the same level. 

Note that the unbound states in such partial definition in the form (42) were introduced 
early (perhaps, for the first time by C. V. Sukumar in [8], see p. 2925-2934) but, mainly, 
when the energy of factorization (when it has real value) locates not higher then lowest energy 
level E^^ of Vi (for example, see [16], [6] p. 323-325 — when Vi has discritc energy spectrum 
completely, see [15] — when Vi has energy spectrum with discrete and continuous parts) or 
has complex value (for example, see [79, 80]). Of cause, an important progress has been 
achieved by B. N. Zakhariev and V. M. Chabanov in development of the main formalism of 
the unbound states: in [49] (see p. 286-290) — for inversion of 5-potentials {£ < 0), in [49] 
(see p. 293-297) — for analysis of deformations and for finding out their rules for rectangular 
well and oscillator after inclusion of one new unbound state (at 8 < £"1^^), in [49] (see p. 309- 
318) — for generalization of the main formalism of such states into a multicannel case, in 
[62] (see p. 378-382, 382-389) and in [81] — for generalization of the formalism of double 
Darboux transformations using unbound states with complex energies of factorization (for non- 
hermitian Schrddinger operators with complex potentials and real energy spectra), allowing to 
shift selected levels in spectrum into complex region of values. In complex consideration of 
solutions for superpotential, potentials (when energy spectra have real or complex values), it 
needs to point out such intriguing direction in SUSY QM as PT-symmetric quantum mechanics 
(for example, see [82, 84, 83, 80, 85, 86, 87, 88, 89, 91, 75]). 



4.2 A general solution of wave function of unbound state 

The second partial solution of wave function of the unbound state for the potential V2 at the 
level can be written down as (it is taken from [6] p. 323-324, as a generalization of a 
solution of WF of the ground state at level E^'^; in contrast to [6] we define the unbound state 
at arbitrary level E^^): 

^^'^^ = (44) 



U^) = J {^^\^')fdx'. (45) 
—00 

Indeed, substituting such function into the Schrodinger equation with the potential (41), we 
obtain equality at the level -E^-*. Let's analyze, what boundary conditions this function satisfies 

to. As (p^^\x) describes the bound state, the function Im{x) is finite at any a; (at a; ^ +00 
one can connect the function Im{x) with a finite normalizing constant for v^^-*). Therefore, at 
a; — > or a; — > a the second partial solution <^^'^'' tends to infinity, i. e. it describes the unbound 
state at E!^\ 
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Let's compose a general solution of WF of the unbound state from two partial solutions (42) 
and (44): 



(2) 



(46) 



Taking into account (46) and (43), we write the superpotential as: 



dx 



dx 



m) 1 



(47) 



/^Wm{Xm) = log \Im + 



Conclusion: 

Changing the parameter A^, one can deform the superpotential Wm-i but it does not change 
the form of the potential V2. Prom here we obtain: 

• WFs of all hound states and arrangement of levels in the spectrum for the potential V2 
are not changed in variations of the parameter Am/ 

• We obtain the whole set of new potentials Vi(Am) (which can be considered as deformations 
of the given V\ by the parameter Xm), which all are partners to one V^. 

Note that the general solutions for WF for unbound states in the formalism similar to (44)- 
(46) (perhaps, which was introduced for the first time by C. V. Sukumar in [8] p. 2926) can 
be found in literature also (one can meet with it essentially more rarely then with the partial 
solutions in the form (42), but which improves essentially the SUSY QM methods of deforma- 
tions) , but again mainly, when the energy of factorization locates not higher then lowest energy 
level E^^^ of V\ (for example, in [6] p. 323-326 — the most universal and convenient variant, 
in [15, 16] — a very improved formalism with analysis of different peculiarities described in 
details, in [49] p. 286-288 — a simple original variant with use of needed boundary conditions 
and its application for inversion of the given potential Vi and in [49] p. 309-313 — a general- 
ization of such formalism for the multichannel case for some special types of potentials, in [48] 
p. 1571-1572 — use of unbound states for insertion of new levels into spectrum in the formalism 
of the inverse problem). 



5 Isospectral potentials 
5.1 1-parametric potentials 

Under changing the parameter Am, the solution (46) keeps interdependence between WF of the 
unbound state for V2 and WF, dependent on Am, of the bound state for Vi at the level -E^-*, 
which coincides with the energy of factorization E. Thus, we obtain a possibility to deform the 
potential V\ and all its WFs (without a displacement of levels in the spectrum) with the help 
of one parameter Am at the selected energy of factorization. The deformed WF at for V\ 
has a form: 

tnl„' Im -r Am 



ifim 
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Let's find the deformed potential Vi{Xm), using its connection witli Wmi^m) and taking into 
account the deformation: 



dWmiXm) 1 



Vo 



Vi + 2 



dW,, 
dx 



dx 



> ^ V,{Xm)=V,-2-^{WU>^m)-Wj)=Vi-2^AWm{Xm), 



and, therefore: 



ViiXJ 
AVi{X, 



V, + AV,{XJ, 

-2^AVF^(A^) = -2^1og 



Im + Xr, 



(49) 



Prom (48) and (49) one can see, how using the parameter to deform WF of the bound 
state at the level E^^ , which coincides with the energy of factorization, and the shape of the 
potential Vi. 

Note, that we have obtained 1-parametric set of new exactly solvable isospectral potentials 
Vi{Xm) on the basis of deformation of wave function of the arbitrary bound state, in contrast 
to [6] (see p. 324-326), where for the deformation the wave function of the ground bound state 
is used only. Such a generalization allows us further in frameworks of the SUSY-approach to 
construct all pictures of deformation of a rectangular well and its WFs (without a displacement 
of levels), which were obtained early by the approach of inverse problem (see [47], p. 916-927). 



5.2 Deformation of wave functions of other bound states 

Now let's find, how the WFs of other bound states for Vi at levels, which do not coincide with 
the energy of factorization S, are deformed by the change of the parameter A„j. WFs of the 
potentials Vi and V2 for coinciding levels E^^ = E^J 7^ ^ = E^^ are connected so (see also [6], 
p. 287-289): 

^(2) ^ IfijCl) _ ^(1)1-1/2^ ^(1) 

(1) ^ |£;(i) _£;(i)|-i/2^+(.(2) l^Uj 

where 

= + ^^^^ 

Here, in calculations of the normahzing constants for these WFs we use a condition of their 
normalization (as for bound states). Taking into account, that WFs of the potential V2 are not 
deformed under variations of the parameter A^, we obtain: 

<HAJ = \E^^}-Ei^)\-^AU\m)Arn^^^} = 





^(1)1- 


= l^iV- 






£;(i)|- 

m 1 




^(1)1- 
m 1 




ni 



dx 
d 



1^^,' - E^::^\-'At,A^^i^) + |E« - Ei^)\-'AW{X^)A^ ^. 



(1) 
ni 



or 



(52) 
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where 

A<)(A„) = ^fi^A^;, 

~ I (53) 
AW{Xm) = ^ log |/™ + A™ . 

Expressions (52)-(53) describe, how under change of the parameter WF of the bound state 
for Vi is deformed at the arbitrary level E^^, which does not coincide with the energy of 
factorization S. 

Analysis: 

• Analyzing a possibility of the potential Vi to change its shape with keeping the shape of 
the potential V2, we come to the following assumption: Each separate level of the energy 

spectrum contains one degree of freedom for deformation of the potential and its spectral 
characteristics (expressed mathematically through the parameter Xm). 

• From (53) we find: The distance between levels E^^ and S = E^} is larger, the influence 
of the parameter \m on the deformation of WF at the level E^^ is smaller. 



5.3 n-parametric potentials 

One can generalize the method of deformation of the potential Vi without displacement of the 
levels in the spectrum by use of the parameter of deformation A^, described in sec. 5.1: in the 
beginning we fulfill a transition from the given potential Vi into a new potential V^. by applying 
SUSY-transformation many times, which transforms consistently the bound states for selected 
n levels for the potential Vi into unbound states for the new potential l^j, and then we fulfill an 
inverse transition from the obtained potential Vn into a new deformed potential Vi (which will 
depend on n parameters Aj, i = 1 . . .n), transforming consistently the found unbound states 
into the bound states (in the assumption, that they exist). By such transformation we obtain 
a whole set of new exactly solvable potentials with identical energy spectra, which depend on 
n parameters of deformation and one can name them as n-parametrical family of isospectral 
potentials. Note, that here we use arbitrary numbers of the levels of the potential Vi and their 
sequence for the n-multiple application of SUSY-transformations (in contrast to [6], p. 326-328, 
where the lowest n levels are used consistently only). 

Let's consider a construction of the family of 2-parametrical isospectral potentials (as it is 
done in [6], p. 326-328 for two lowest levels). Let's we have a potential Vi with known WFs 
and energy spectrum. In such spectrum we choose two levels E'^j and E^^ with numbers mi 
and In the beginning we fulfill a transition from the potential Vi into a new potential 
V2 applying SUSY-transformation, for which we define a superpotential W^l with energy of 
factorization Srm, which coincides with (as it is done in sec. 5.1). As we know, at this level 
there is WF of the unbound state for the new potential V2. 

^^KAi) = %y^, 1 {^t\{x')fdx\ (54) 

99^^^ is WF for Vi at the level and Ai is an arbitrary constant. One can find WFs of the 
bound states for other levels of the potential V2 from the first expression in (50) and from (51) 
at m = mi, where: 

^l^l - (55) 
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and one can obtain a shape of the potential V2 from: 



V2 = V, + 2^iy(\) = V^i - 2^ log |. (56) 

Further, we fulfill the transition from the potential V2 into a new potential V3 applying next 
SUSY-transformation, for which we define the second superpotential VF^] with the energy of 
factorization £,^2^ which coincides with another level E^^. For the new potential V3, we obtain 
WF of unbound state for this level E2,1: 



TO2 ■ 



= (57) 

and WFs of bound states for other levels: 

^i? - \E^^^ - E^^^-'/'A^^l^^S at E^) = (58) 

where 

J{^mx')fdx', Agl^l + W^^l P^(^i = -^log|^(^)| (59) 
—00 

and a new constant A2 is introduced. We calculate the shape of the potential V3 as: 
One can see, that a function 

describes unbound state for the potential V3 at the level E^^^ also (here we have introduced a 
"normahzing" factor). We see, that consecutive n-multiple application of SUSY-transformation 
to the given potential Vi at the selected n levels reduces number of the hound states by n and 
increases number of unbound states by n also. 

Now wc fulfill an inverse transition, coming back from the potential V3 to a new potential 
Vi with consecutive transformation of the unbound states into bound states at the levels E^_^ 
and However, in definition of superpotentials for two inverse SUSY-transformation we 

can change sequence of a choice of the energies of factorization and £^2 (in contrast to 
[6], p. 326-328). Let's introduce new designations fhi and Aj for numbers of the levels and the 
parameters by such condition: 

(mi = mi,m2 = m2, Ai = Ai, A2 = A2) or (rrii = m2, m2 = mi, Ai = A2, A2 = Ai). (62) 

For definition of the energy of factorization (let's denote it as in the construction of the 
superpotential for the transition from V3 into V2 we choose the level E^^ with the number mi. 
The unbound state (/9^j|(Ai) of the potential V3 at this level transforms into bound state v'^|(Ai) 
of the potential V2, dependent on one parameter Ai (a question of existence of divergences is 
important in obtaining of the solution for WF of such "bound" state; but, perhaps, one can 
exclude them as it will be done in sec. 6 in deformation of a rectangular well with infinitely high 
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external walls with the energy of factorization, located not below than the first exited level of 
the spectrum of such potential): 



rm2 



^ml (A2) 



J-m2 ~r A2 



The other unbound state </^^2('^2) of the potential at the level E^^ with the number m2 in 
the transition forms a new unbound state (p^^(Xi, A2) for the potential V2: 

For WFs of bound states for the other levels we obtain: 
where 



at mi = m2, 
at fhi — mi. 

(2) 



(63) 



I (2), 



d 



d 



.(2), 



d_ 

dx 



log|<^S(Ai 



The new potential V2(Ai) already depends on one parameter Ai and has a form: 



y2(Ai) ^V,- 2^W^4?(Ai) = ^3 + 2^ log |^g!(Ai) 



(64) 
(65) 
(66) 

(67) 



For the last transition from V2 to Vi , for the definition of the energy of factorization Sf^^ for 
the second supcrpotential we choose the level with the number 7712. In result, we obtain: 

• the unbound state of the potential V2 at the level E^l transforms into bound state of the 



potential Vi 



^PrnliXi, X2) 



1 



<^S(Ai, A2) 



(68) 



• the bound states of the potential V2 at the other levels transform into the bound states 
of the potential 14 : 



<nAi, A2) = \EW - Egl\-'/'Agl^{X,, X^W^i), at E^ = Eg 



112 ' 



(69) 



where 



>lS^(Ai, A2) ^-^ + W^S(Ai, A2), wl'i(Ai, A2) = log \ipglCXi, A2)|; (70) 



d 



the new potential Vi already dependes on two parameters Ai and A2, and has a form: 



K(Ai,A2) = y2(A2) -2^PF4'i(Ai,A2) = l^2(A2) + 2^1og|(^S(Ai,A2)|. (71) 



Taking into account (56), (60), (67) and (71), we obtain: 



T4(Ai,A2) =T4-2^1og 



^(1)^(2) 
T nil ' ff^2 



<^mi(Ai)<^r;^2(Ai, A2) 



(1), 



(72) 
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Let's generalize the found result into a case of deformation of the potential Vi by use of n 
parameters Aj. For this purpose, in the beginning we choose n levels in the energy spectrum of 
this potential. Let's these levels are numbered by two independent sequences of numbers mi, 
m2 . . . rrin and mi, 7712 ■ ■ ■ m„, accordingly. Then for the deformed potential we find: 



^i(Ai...A„) = yi-2^1og 



(1) 



(-^1)95)^2 '^(-^1, A2) . . . (f)^i{Xi . . . A„) 



(n-l) 



(1) 



or 



Y^mi ■ ■ ■ T'lUn rmi 



(Ai)(^g(Ai,A2)...<^^i(Ai...A„) 



(2) 



(73) 



(74) 



WFs . . . ifl^lj^ of the bound states for the potentials Vi ■ ■ - Vn can be found on the basis of 
. . . i^')^}^ of the bound states of the started potential V\ (we write new WFs up to normalizing 
constants, denoting them as C): 



^ mi r m2 ' ' 
C^(2)^(l) (1) 



mi 
(n) 



(75) 



The deformed WFs </7^. ^"^""^^ ( Ai . . . Aj) of the bound states with the parameters Aj can be cal- 
culated on the basis of WFs (p!f^^^^ ■ ■ ■ <^^^^^ of the unbound states of the potential Vn+i (at 
i = 2...n): 

VrnAXl)- , 



(n-l) 



V-'mi 

(Ai,A2) = 



(n-i+l) 
rrrii 



f^^niXi, A2 . . . A„) 



(Ai, A2 . . . Aj) 



1 



(76) 



(n-i+2) 



(Ai, A2 . . . Aj) 



-(2) 



or 



¥'S(Ai, A2) 



(Ai, A2 . . . A„) 

CA^}'^{Xi) ^m^^\X2 



'i^in-.i (^li^2jA3) 

(Ai . . . A,) 



{n-i+2) 
rrii 

.(2) 



cAi;'>'^x^,x,)4i^{x,) ^i;:r>{Xs) 



,(n-'i+2),+ . 



^fflniXl . . . A„) 



CA)^._'^^''^{Xi, A2 . . . Aj_i) 
C^kl'^L (Ai, A2 . . . A„_i) 



• • ^mi' (Al) 

^S'+(Ai) (^kr^(A„ 



(77) 



WFs of the unbound states of the potential Vn+i have forms: 



(n+l) 
m2 



(A2 



= CA(;:l...A(j^0gl{x,), 



(78) 



WFs (/'^■^^^ (Aj) must be calculated by such formulas (at i = 1 

+ A 



(i) 



n): 

(v.»(x'))'cix'. 



(79) 
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The non-deformed and deformed operators A and A'^ have forms (at i = 1 . . . n): 

= ^ + ^-^^"'"""(^^ • • • = + Wkr^'\x, . . . AO, (80) 

where 

W^] = log I, l^4r+^)(Ai • • • A.) = log l^kr^^^Ai . . . A.)|. (81) 

The deformed WF of arbitrary bound state for the potential T4(Ai . . . A„) at the level with 
the number where ^ E^^^, . . . E^ ^ E^^, has a form: 

^«(Ai ...K)= CAgl:^{X, ...K)... A^l'^CX,) . . . . (82) 



5.4 New approaches of construction of isospectral potentials 

It turns out, that if in definition of the function of factorization for construction of the su- 
perpotential to do not restrict oneself only to the bound states of the given potential Vi, but 
to use an arbitrary unbound state with a selected level for this potential, then such a way 
allows to construct new types of potentials, which can be isospectral to Vi. Here, use of the 
unbound states is caused by rejection of imposing boundary conditions (for bound states) on 
WF of such state, that does not destroy interdependence (but, rather, expands it) between 
two potentials-partners. Such consideration inevitably leads to an absolute arbitrariness in a 
choice of a value of the energy of factorization ( which can be not coincident with the levels of 
the energy spectrum of the given potential Vi). Let's consider this question in more details. 

Let's we know the wave function (p^^ of the unbound state at the level E^^ for the given 
potential Vi. Then a function of a form: 

/r(AJ = ^SL±^, /(i)(^)= l[^^^\y)fdy (83) 

is a general solution of the Schrodinger equation at the same level E^ for a potential Vq with a 
previous number (which can be named as the number "down"; see [60], p. 448-451), which 
has a form: ^ ^ 

Vo^Vr- 2iy W ^ (^iy(i)^ _ {wi^)^' + E^, (84) 

where 

W^^^^ = ^log|^«| (85) 

and A^ is arbitrary constant. 



Proof: 

For the function 1^^ {x) we have: 



Then we obtain: 




(1)^" 



-(1) I -(1) 1 -(1) I -(1) 



23 



Let's analyze, whether the function fw\^w) satisfies to the Schrodinger equation with the potential Vq at the 
level Eyj. We write: 

Taking into account (86), we obtain: 

One can see, that this expression becomes identity, if to define the potential Vq as (84) with taking into account 
(85). Therefore, we have proved, that the function fw^ in the form (83) is the solution of the Schrodinger 
equation for the potential Vq at the level E^j- One can see also, that the solution (83) consists of two linearly 
independent partial solutions and, therefore, it is a general solution. 



Consequences: 

• The potentials Vq and Vi are SUSY-partners with the superpotential of the form (85). 

• Variations of the parameter with the given form of the wave function (p!^^ do not 
deform the superpotential and, therefore, the potential Vq. 

• Prom (83) one can see, that by use of the parameter one can deform the function 
(with keeping of the shape of the potential Vq). 

• The found solution for f^^ describes the bound state of the potential Vq at the level Eyj 
not always, but only in the case, when it satisfies to the boundary conditions for the 
bound states and a condition of a continuity of this function inside whole region of its 
definition (which can be found from an analysis of possible divergences of this function). 



In dependence on whether such a value of the parameter exists, at which the function 
/^°^ satisfies to the bounding conditions, we have two cases (if the energy of factorization Eyj 
does not coincide with any level E^^ of the spectrum for Vi): 

• If the function f^^ satisfies to the bounding conditions, then it defines a new bound state 
for the potential Vq at the level Eyj. Then one can consider this function as WF of the 
bound state for Vq (let's denote it as </^^^), and the spectrum of the new potential Vq 
differs on the spectrum of the given potential Vi by the additional new level Eyj. 

• If the function f^^ does not satisfy to the bounding conditions, then it defines a new un- 
bound state for the potential Vq at the level E^ (let's consider this function as WF 
of the unbound state and denote it as <^2i''')- Then if the superpotential for SUSY- 
transformations keeps all other states for Vq as bound ones, then the spectra of the 
potentials Vq and Vi are equal absolutely, and these potentials are isospectral. 



24 



The second case gives us a new way of construction of the isospectral potentials on the basis of 
the known WF of the unbound state of the potential Vi (note, that this way is some easier in 
comparison with the approach from sec. 5.1, however here it needs to pay more attention to an 
analysis of divergences). In the new approach one can deform the found isospectral potentials 
Vq not only by variation of the parameter A^,, (as in the approach from sec. 5.1), but also by 
change of the value of the energy of factorization in a possible range and by deformation of 
WF shape of the unbound state (^^^^ . 

Now let's redefine the superpotential on the basis of the function f^\Xw) so: 



By use of this superpotential one can construct the new potentials Vo(A„,) and Vi(A„,). These 
potentials will be SUSY-partners and, therefore, they should have the energy spectra, which 
coincide completely with the energy spectra of the old potentials Vq and Vi (with a possible 
exception of the level i?^). 

The definition (87) introduces the dependence on the parameter A^ into the superpotential 
(and, as a result, into the potentials Vo(A^) and Vi(A^)). In result, we obtain the second 
new approach of the construction of new isospectral potentials, which can be some like to the 
method from sec. 5.1. A difference between these two methods consists in the following: 

• If the method in sec. 5.1 is constructed with use of the superpotential, which is defined 
on the basis of WF of the bound state of the starting potential Vi, here we define the 
superpotential on the basis of WF of the unbound state of the starting potential Vi. Such 
a way is richer essentially and it gives more large class of new isospectral potentials. 

• In contrast to the method in sec. 5.1 for the construction of the new isospectral potentials 
with the transition of their numbers "up-down (Vi ^ V2 Vi(A)), here we fulfill the 
inverse transition of the numbers "down-up (Vi Vq ^ ^i(A^)). 

Conclusions: We have obtained two new methods of the construction of the set of the isospec- 
tral potentials on the basis of the given WF of the unbound state for the given potential V\ 
at the given level E^. In contrast to the known approach in sec. 5.1, where as the parameter 
of deformation we use \, in the new approaches as two additional parameters of deformation 
the function of factorization (defined on the basis of WF of the unbound state) and the energy 
of factorization E^, are used, that makes a class of the new isospectral potentials larger, and 
methods richer. 

The first approach: 



d 

dx 



log|/r(A^)|. 



(87) 




(88) 
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The second approach: 



AV^{E^, ^W, A^) = -2^ log | + A^|, 

X 

— CX3 

— the parameters of deformation. 



(89) 



Consequence: On the basis of the second approach one can generahze the method from sec. 5.3 
for the construction of the n-parametrical family of the isospectral potentials with taking into 
account of a possibihty to define the function of factorization as WF of the unbound state. 

6 A rectangular well with infinitely high walls 

An important progress has been achieved in development of methods in the approach of inverse 
problem for construction of new exactly solvable potentials on the basis of known ones. One 
can say, that today the methods of the inverse problem come in the vanguard in forming an 
unified theory of exactly solvable models (sec [47], p. 916). Thus, in [47, 49] a nice analysis of a 
deformation of a rectangular well (and also some other potentials) under variations of spectral 
characteristics, allowing to observe an influence of these characteristics on the potential shape, 
was fulfilled on the basis of such methods. One can explain an interest to such papers by words 
of authors (see [47], p. 916): "Approach of inverse problem is remarkable by that it allows from 
another side to see on connection of forces with physical characteristics of systems, where they 
act. A possibility appear to construct the systems with needed properties..." (in the original 
text: " , , . ...")• 

It turns out, that supersymmetric methods, developed at full their power, allow to achieve 
this also, not yielding to the methods of the inverse problem and having own attractive aspects. 
Let's consider, how one can deform the rectangular well with infinitely high walls (which we shall 
consider as a starting potential V^) by use of the method from sec. 5.1, where in construction 
of superpotential an energy of factorization coincides with arbitrary selected level of discrete 
spectrum of the given Vi, and a function of factorization is FW of a bound state at such level. 

Let's define this potential so: 

1 0, atO<^<a; 

[ +00, at X < U or X > a. 

The problem of determination of the energy spectrum and WFs for this potential is, perhaps, 
the task, which is included into any textbook at quantum mechanics. Solving Schrodinger 
equation with the potential (90) at the energy E higher then a bottom of the well {E > 0), we 
find a general solution of this equation: 

(p^^^ (x) = A sm{kx) + B cos{kx) , (91) 

where A, B are arbitrary (complex) constants and k = \/E. 



26 



For extracting from the solutions (91) WF, which describes the bound states, we use bound- 
ary conditions (providing continuity of WF in the whole region of its definition) : 

(^W(0) = and (^W(a) = 0. (92) 

Application of the first condition gives B — 0. From the second condition we obtain: 

{ka ^nn)^[k = kn^^^ ^m) ^ ^ kl = ^n' = E^n''^ . (93) 

at n = 0, 1, 2 . . .. So, we obtain the discrete levels of the energy spectrum E and then to each 
level with its WF we add corresponding index n. The constant A for WF with the number n 
can be found from the condition of its normalization: 

^^:\x)\'dx = 1, (94) 



From here we obtain: 



J \(p^^\x)\'^dx = J \A\'^\sm{knx)\'^dx = lAl"^ J -{1 - cos{2knx)fdx ^ 







\A\' _ sin{2knx) 



2kn 



x=a 



x=Q 



\A 



2 



-a = 1 



or 

'a (^^' 

We see, that the coefficient A is determined up to phase factor and we choose: A — \A\. This 
gives us only real values for WFs for all states. Note, that the coefficient A is the same for 
different states with the numbers n. Finally, write the solution for the normalizable WF of the 
bound state with the number n: 

ip^\x) = \ -sm{knx) = \ -sm{kinx). (96) 

yd yd 

From (93) and (96) we find coordinates of nodes for WF with the number n: 

xi^^^-, 1^0,1,2. ..n. (97) 

Kn ri 

One can see, that the nodes for WF with arbitrary number n are located on equal distances one 
from another, which decrease in increasing of the number n. 

6.1 The deformation of the rectangular well 

Let's find a function I^{x): 

4'^^)= / {^^\x'))Hx' ^ - j sm\krnx)dx' ^ - j -[l - cos{2k^x'))dx' ^ 

-oo 



[ , _ s\n{2kmx') 



a V 2kr, 



^ ^ 1 / siii{2kjnx) 

x'=0 ^ '^^m 
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or 



CI \ Qikfji 



(98) 



For the deformed potential T4(A^) we obtain: 

Vi{x, \m) = log \I^ + A^l = log 



-\x 
a 



sm2kjnX\ 



or 



Vi(x, A^) = -2— log 



(99) 



Prom here one can see, that at A^ — > +oo (or A^ >> 1) the potential Vi{Xm) tends to the 
starting undeformed potential Vi. 

We find an explicit form of the deformed potential Vi{Xjn)- Prom (99) we write: 



Vi{x,Xm) = -32fe 



sin kinX (^A:,„(,T + aA„j) cos k,nX — sin k^x^ 
{2k jYiX sin 2kfjiX + 2kjnO'Xf}ij 



(100) 



Le^'s name points, where the potential Vi equals to zero, as zero-points or "nodes" of this 
potential. Then according to (100), all zero-points of the deformed potential Vi(Am) can be 
separated into two types (for the rectangular well such a property has found for the first time) : 

• Zero-points, which coincide with the nodes of WP of the bound state at the level with the 
number m, coincident with the energy of factorization. These zero-points do not displaced 
under deformation of the potential by use of the parameter A^ and their coordinates equal 
(they determined from (96) with change of the index n into m): 



ttI al 

= — , / = 0,l,2...m. 



k. 



m 



;ioi) 



Zero-points, which do not coincide with the nodes of WF of the bound state at the level 
with the number m, coincident with the energy of factorization. These zero-points are 
displaced under the deformation of the potential by use of the parameter A^ and their 
coordinates are determined from: 



(102) 



6.2 The superpotential under the deformation 



Let's find the superpotential for the given potential Vi. 



d I fill dsinkmX 



Sill (^kfYi^j^ — ^Tfi cot JZfjiX 



or 



Wm{,x) = —krnCOtkmX. (103) 

Note, that the superpotential has divergences at the nodes of WP with the number m (see (96)). 



28 



6.3 The deformation of wave function of the bound state at the 
level, coincident with the energy of factorization 

For the deformed potential Vi(Am) we find WF of the bound state at the level with the number 
m, coincident with the energy of factorization. The expression (48) determines such function 
up to a normalizing constant. One can write: 



d kfjiX sm IvfjiX cos k^yiX ~\~ k')yi(i\fyi Jv^iX siii Jv^yiX COS k^p^x ~\~ kffiQj\.jyi 

or 



^p{^)(^x X ) - C \/2a/c^sin(fc^a;) 

kjjiX sm kjYiX COS kjyiX ~\~ kfjiCi\ui 



(104) 



or 



The constant can be found from the normalizing condition of the found WF: 

a 

I |(/p«(x,Aj|2dx = l (105) 



1^ |2 2ak'^sm'^{kmx) _ 

( kmX sin kjjiX cos kjjiX + k^nClXxn j 

\Cm\'^ ■ 2(ifc^ j — '- ^dx 

1^ |2_ ,2 f 1 - COS (2kmX) 



Substitution: 



y — ^m,3j 2 sin 2A;mX + /c^oA^, (x — 0) > y — k^aX^^ 

dy = km{l - cos 2kmx)dx, {x = a) {y = kma{l + A^)). 



Then we obtain: 



x=a , 

' x=a /I 2 



\r l'^ . nk I ^ — -\r \^ . nk - 

\ y y 
1 1 

I C*™, I ( ; : — T ] = \C, 



x=Q xkjnttXjTi kmO'iXm ~\~ 1) 



-Xm A^ + iy Aj„(A^ + 1) 

From here we find: 

\Cm\ = Va^(A^ + 1). (106) 

We see, that the normalizing constants Cm are different for the different m and A^. From (104) 
and taking into account (106), we obtain a final expression for the normalizable WF (we select 
Cm — \Cm\)'- 



^(1)^^ A ) — ^Jx (A + 1) v2afem sin(A:m.a^) (107) 

kfjiX sin kf-f-fX cos kjyiX ~\~ Aj^ciA^ 
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Property 1: nodes of this WF do not displaced at its deformation by the parameter Xm for 
arbitrary m and they coincide with the nodes (97) of the undeformed WF. Let's find a derivative 
of this WF: 



— \J Xm{Xm + 1) 



dtx dx JZfjiX sm kuiX cos kfj^x ~\~ ^ttt.gE'A^ 

/ COS kvyrX ( kjj-iX sm kynX cos kjj-,x ~\~ kjjiCi\'f-n } sm kyjiX i 1 cos kjj-,x ~\~ sm kynX ) 

- ^A„(A^ + l)V^kl ^ ■ 1 ^-2 ^ 

( kfYiX sm kfjiX cos k^iX ~\~ A^tt^^^A^ ) 

nr-. — 77 77, o cos kfYix(^kfjiX sm kfYfX cos kuiX ~\~ kfjid^ui^ 2 sm kuiX 

— Y ^ay^mi^m + J-j/^m 7 \ ^"2 

( kfYiX sm kjjiX cos kyyiX ~\~ kyyiQj\f- 



or 



d^rni-^t ^rn) 7 77 ' 77, 2 ^'^^ kmX{kjnX — Sin kjnX COS + kmdXm) ~ 2 sin"^ kmX 

— y ZaAm['\m + J-jfe 



y -Illy -III ' /■-m / \ z 

( kffiX sm kfYiX cos k^^x ~\~ kfuCiXfYi | 

(108) 

At x-coordinates x — and x — a of the walls we obtain (cos/c^a = (—1)™): 



d^^K^Am) ^ / 2(1 + A^) c;99W(q^aJ ^ {-l)'^km / ^^"^ (109) 
dx '"V oA^ ' dx ^\a{l-\-Xm) 

Property 2: According to (109), angles of slope of WF relatively axis x at x-coordinates of the 
walls are changed at the deformation of this WF by use of the parameter Xm- 

6.4 The deformation of wave functions of other bound states 

According to sec. 5.2, wc write the normalizable WF of the arbitrary bound state for the 
deformed potential Vi{x, Xm) at the level with the number n, not coincident with the energy of 
factorization {En 7^ ^ = Em) '■ 



^Vn \^m) — 7— 7n ^7777^' 



,^^^-s(!^r"^'^"^' (110) 

^W{Xm) = ^log|/^ + A^|. 

We see, that the deformation of this WF in result of change of the parameter A^ is defined by 
the function Aip^^^(Xm) purely. Taking into account (98), (103) and (96), we find: 

AW^(AJA^(^« = log \Im + Xm\(-^ + Wm)v^^rl^ = 

cix \OjX / 

4A:^sinA:^x / , . , , i, ■ i, \ 

; ~z ~z ~ 1 Kfi COS fbijX sm Km^ '^jxi COS Km-^ Sm K^X J . 

CL Zfvm-^ — sin Zfvm-^ ~t~ ^km^Xm 

From here it follows, that the function Acp^^Xm) can be equal to zero only at fulfillment of 
one from two following conditions: 

sin kmX = 0, kn cos knX sin kmX — km cos kmX sin knX. (112) 
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The first condition is fulfilled at such points: 

xi = ^ = —, at Z = 0,1,2... m. (113) 

Thus, we have obtained the following properties of the deformation of WF of the bound 
state at the level with arbitrary number n, not coincident with the energy of factorization 
{n 7^ m, for the rectangular well these properties have obtained at the first time). 

Properties 1: 

• There are such points, where all deformed curves of WF of the state with the selected 
number n intersect between themselves and with undeformed WF. One can separate these 
points into two types: 

— The points of the first type, x-coordinates of which equal to (113) and coincide with 
the nodes (97) or (108) for WF of the bound state at the level with the number m, 
coincident with the energy of factorization, and with zero-points of the potential of 
the first type (101) also. 

— The points of the second type, x-coordinates of which do not coincide with the nodes 
of WF of the bound state at the level with the number m, coincident with the energy 
of factorization. The coordinates of these points can be determined from the second 
equation of the system (112). 

• The boundary points a; = and x — a are such points of both the first type, and the 
second one. 

• Such point of the first and second types do not displaced at the deformation of WF using 
the parameter A^. 

• WF of arbitrary bound state (in its deformation using the parameter A^) have m+1 points 
of the first type and number of the points of the second type, which can be determined 
from the second equation from (112). 

Now let's find a derivative of this WF. From (110) we obtain: 

d^n\>^m) ^ ^ — A^«(A^) = 

dx dx dx ^ ^ 

difffl^ /2"4/c^ cos kmX (kn cos knX sin kmX — km COS kmX sin knx) 

dx \l a I Ei^^ -E^m^l (2kmX - sin 2kmX + 2k^a\m) 
2 4:km (km ~ ^n) si^^ siu knX f2 16/of„ sin^ kmx[kn COS knX sin kmX — km COS kmX sin knX 



a \E^^^ - E^m^ I (2kmX - sin 2kmX + 2A;^aA^) V a |^a) _ ^(i) | (^j^^^ _ ^^^^ ^kmX + 2A;„aA^)^ 

(114) 

From here we have in the coordinates x = and x = a oi the walls: 

M;(M = Mi, i!-Arf)(A„) = 0. (115) 

dx dx dx 

Thus, we have proved the properties, pointed out in [47] (p. 918): 
Properties 2: 
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• At the coordinates x = and x = a oi the walls WF of the arbitrary bound state at 
the level with the number n, which does not coincide with the energy of factorization 
{En ^ — Em), leaves at the same angle of the slope to the axis x under its deformation 
using the parameter X^, (or under change of a derivative 7^ in the approach of the inverse 
problem) . 

• The angles of leaving of this WF from x-coordinates of the walls are determined by the 
derivative of the undeformed WF in the form (115). 



6.5 Divergences in solutions 

The found solution (108) for WF at the level with the number m, which coincides with the 

energy of factorization, or the solution (110) for WF of the state with the number n at another 
level for the deformed potential describes the bound state only in such a case, when this WF 
has not divergences and is continuous on the whole region of its definition at x. According to 
these solutions, the divergences for WF with the selected numbers m and n appear only in the 
case, when: 

2kjnX — sin 2kjnX + 2kmCiXm — 0. (116) 

One can exclude these divergences (at any x) from the found solutions for WF, in particular, 
by imposing on a choice of the parameter the following requirements: 



or 

^^^_l + 2Kna^ A„>-^. (117) 



2 A/777, CI 2 A/777 ot 



From here we conclude the following. 
Properties: 

• If the parameter is in one from regions, determined in (117), then WF of the form 
(108) for the state at the level with arbitrary number m, which coincides with the energy 
of factorization, does not have the divergences, it is the continuous function at < x < a 
and equals to zero at coordinates x — and x — a oi the walls. Therefore, it describes the 
bound state with the number m (in spite of the fact that the superpotential, constructed 
on the basis of WF of the excited bound state with the number m, has the divergences 
at the nodes xi — al/m of this WF, i = 0, 1, 2 . . .m). 

• If the parameter A^ is in one from the regions, determined in (117), then WF of the form 

(llO)-(lll) for the state at arbitrary level with the number n, which does not coincide 
with the energy of factorization, docs not have the divergences, it is the continuous 
function at < a; < a and equals to zero at the coordinates a; = and a; = a of the 
walls. Therefore, it describes the bound state with the number n (at arbitrary m, in 
spite of the fact that the superpotential, constructed on the basis of WF of the excited 
bound state with the number m, has the divergences at the nodes xi — al/m oi this WF, 
I = 0,l,2...m). 
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• If the parameter is in one from the regions, determined in (117), then the deformed 
potential of the form (100) does not have the divergences and is the continuous function 
at < X < a (in spite of the fact that the superpotential, constructed on the basis of 
WF of the excited bound states with the number m, has the divergences at the nodes 
xi = al/m of this WF, Z = 0, 1, 2 . . . m). 

Conclusion: 

We have proved an appropriateness of the method (from sec. 5.1) of construction of new isospec- 
tral potentials (with own hound states and without the divergences) on the basis of one known 
potential Vi in the form of the rectangular well with infinitely high external walls, if in the 
definition of the superpotential as the energy of factorization the arbitrary selected level of the 
discrete energy spectrum of the given V\ is used, and as the function of factorization the WF 
of the bound state at such level is used. 

6.6 Analysis 

Let's consider, how in the frameworks of the supersymmetric approach one can deform the rect- 
angular well with infinitely high external walls and finite width a (further, for the convenience 
we choose a — 1). As a method of deformation we shall use the method from sec. 5.1, where in 
construction of the superpotential as the energy of factorization we choose an arbitrary selected 
level with number m of the energy spectrum of this well, and as the function of factorization 
we use WF of the bound state at this level. We shall consider each new potential, obtained by 
help of this method, as a deformation of the old rectangular well (i. e. the starting potential 
Vi{x, \m) in its formalism) using parameters of the deformation — the number m of the level, 
equals to the energy of factorization, and the parameter A^. We shall obtain the deformed 
potential from (99), the superpotential — from (103), WFs — from (108) and (110). 

The found deformed potential (a) and its FW of the bound state at the level with the number 
m (b), coincident with the energy of factorization, are shown in Fig. 1 at choice of the different 
values m. In the first figure (a) one can see, that the selected level m determines number 
of intersections of the potential with axis x (or number of "zero-points" of this potential); 
increasing it, one can increase the number of such points of intersections. We also sec, that 
by increasing m both an angle of a slope of the potential in its leaving from the left boundary 
at point a; = 0, and an angle of a slope of WF at its leaving from this boundary, increase 
(that follows from the choice of the definition of the function Im{x) for calculation of the 
deformed WFs and the potential). From these figures one can see, that at the selected values 
of the parameter A^ there are no any divergence on the whole region x G [0, a] both in the 
deformed potential, and in WF at the different values m, and also that WF equals to zero at the 
boundary points x = and x = a (that points out the bounding of this state). It demonstrates 
an appropriateness of the method from sec. 5.1 for construction of new isospectral potentials 
using the superpotential, defined on the basis of WF of not only the ground state atm — 1, but 
also the arbitrary excited state at m ^ 1. 

The deformation of the potential at the change of the parameter Am is shown in Fig. 2. In 
the first figure (a) we see, how the potential is deformed, if as the energy of factorization the 
lowest level at m = 1 is used. Here, we note the following: 

• We see, that this figure qualitatively looks Hkc Fig. 1 (a) from [47] (sec p. 917), in which the 
deformation of such potential by help of the methods of the inverse problem is shown (here, 
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Fig. 1: The deformation of the shapes of the potential and WF of the bound state at the level with 
the number m, coincident with the energy of factorization, in dependence on a change of the number 
m (values m = 1,2,3,4 at Am = 1 are used): (a) — the deformation of the potential; (b) — the 
deformation of WF. 




Fig. 2: The deformation of the potential at the change of the parameter Am (the values Am = 0.5, 1, 2, 5 
are used): (a) — at m = 1, (b) — at m = 2, (c) — the shape of the potential is deformed continuously 
at the change of the parameter Am (without the divergences) 
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varying the parameter Ai and changing the definition of the function /^(x), one can obtain a 
full similarity between these figures). 

• We see, that the increase of the slope of leaving of the potential in the left wall (and the increase 
of minimums and maximums of the potential) at increasing of the parameter Ai coincides with 
the increase of the slope of leaving of the potential in the right wall (and with the increase of its 
minimums and maximums) in Fig. 1 (a) in [47] (sec p. 917), which was caused by a change of 
the derivative 71 of WF of the ground state in the right wall at a; = a. It points out a similarity 
between the parameter Ai in the considered method in the frameworks of SUSY QM and the 
coefficient 71 in the approach of the inverse problem. 

• The following property is fulfilled: the deformed curves of the potential Vi(Am) at the different 
values Xm intersect the axis x at different points, i. e. we obviously observe a displacement of the 
central zero-point of the potential at its deformation by X„i . It is explained clearly (and at first 
time) by the found property from sec. 6.1, and this property is observed also in Fig. 1 (a) in [47] 
(see p. 917) (the considered node is the node of the second type, its coordinate is determined 
by equation (102) at the given Xm)- 

In the second figure (b) we show the deformation of the well, if in the definition of the energy 
of factorization the level of the first excited state at m = 2 is used. Here, we note the following: 

• The deformed well after its mirror refiection coincides practically with Fig. 1 (5) from [47] 
(see p. 917), where the deformation of such well, obtained by the methods of the inverse problem 
at the change of the derivative 72 of WF of the first excited state at the right wall at x = a, is 
shown. 

• One can see, that the increase of the leaving slope of the potential in the left wall, caused by 

increasing the parameter A2, coincides with the increase of the leaving slope of the potential in 
the right wall in Fig. 1 (5) in [47] (see p. 917), caused by a change of the derivative 72. It proves 
the conclusion about that the parameters Xm with the different m in the method from 5.1 play 
the same role in the deformation of the potential, as the coefficients 'jm ( defined at the same 
levels with the numbers m) in the approach from [47]. 

• Prom the figure one can see, that the well after its deformation has "zero-points" of both types, 
described in sec. 6.1 (this property has observed and explained at the first time): 

— three zero-points, located at points of the walls of the well and at its center, which are not 
displaced under the change of A^ (they coincide with all nodes of WF of the first excited 
state at m = 2; coordinates of these zero-points can be found from (101)); 

— two zero-points, located between the previous three zero-points, which are displaced under 
the change of A^ (their coordinates can be found from (102) at the given An^). 

As we see, this property is displayed also in Fig. 1 (5) in [47] (see p. 917). 

According to the figures (a,b), the m is larger, the number of zero-points on the region x G [0, a] 
is larger and the displacement of zero-points of the second type in the deformation of the well 
is smaller. 

Now let's analyze, how the shape of WF is deformed at the level with the number m, 
coincident with the energy of factorization. In Fig. 3 such deformation of WF at the change of 
the number m of such level is shown. From the figures one can see the following: 
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Fig. 3: The deformation of the shape of WF at the level with the number m, coincident with the 
energy of factorization, caused by the change of the number m (the values m = 1, 2, 3, 4 are used): (a) 
— at = 2; (b) — at = 0.1 



• Distances between arbitrary two neighboring nodes of WF with the arbitrary number m are 
equal, they decrease at increasing of the number m and their coordinates can be found from 
(97) (with substitution of index n into m). 

• For any m at decreasing of the module of A^ a relative weight (amplitude) of WF near one wall 
is increased and near another wall is decreased, at Am — > ±cxd a smoothing of the deformation 
of WF takes place and WF tends to its undeformed form. 

In Fig. 4 a dependence of the shape of WF at the level with the number m, coincident with 
the energy of factorization, on the parameter A^ is shown. In the first figure (a) the deformation 




Fig. 4: The deformation of the shape of WF of the bound state at the level with the number m at 
changing of the parameter Am (at Am = 1, 2, 3, 4): (a) — for the ground state at m = 1; (b) — for the 
state at m = 5; (c) — the shape of WF is deformed continuously (without divergences) at change of 
the parameter Am 
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of WF of the ground state at m = 1 is shown. We see, that this figure after mirror reflection 
looks hke Fig. 1 (2) in [47] (see p. 917), where the deformation of WF of the ground state, caused 
by the methods of the inverse problem at change of the derivative 71 of this WF in the right 
wall, is shown! Here, wc reproduce explicitly a demonstration of a deforming property of this 
WF, which was described in [47] (see p. 917-918, and obtained by the approach of the inverse 
problem): "... acceleration of increase o/$i with moving off' from x = a, must be compensated 
by quicker passing to falling, because the integral from the square of wave function of the ground 
state remains normalizable to one. Just the hollow of the potential near x — a provides such a 
behavior ^i{x) ... But, further, it needs to compel the function, with not changing a sign, — 
nodes for the ground state are unallowed, — to tend to zero at point a; = 0. This is achieved 
by reflected part of addition to the initial potential in the left part of the well: under the barrier 
the function is smaller essentially then unperturbed one." (in original text: "... $1 x — a 
, . x^a ^i{x) ... , , — , — X = 0. ; ."). Note, that 

increase of the slope of WF in the left wall (and, correspondingly, decrease in the right wall) 
at decreasing of the parameter Ai in Fig. 4 (a) coincides with increase of the slope of WF in 
the right wall (and, correspondingly, with decrease in the left wall) at changing of 71 in Fig. 1 
(2) in [47] (a possibility to change the angle of the slope of WF at ^-coordinates of the walls is 
proved analytically in sec. 6.3). In the second figure (b) the deformation of the shape of WF 
at the level with the number m = 5 is shown. Here, we note the following: 

• We see, that the deformation of this WF looks qualitatively like the deformation of WF of the 

first excited state in Fig . 1 (7) from [47], caused by change of the derivative 72 in the right wall 
(with taking into account of decrease of m). From here one can see, that the parameters Xm in 
the method from 5.1 play the same role in the deformation of WF as coefficients 7^ (defined at 
the same levels) in the approach from [47] . 

• All deformed curves of WF with the selected number m at the different values of the parameter 

intersect through the same points — the nodes of the undeformed WF of the form (97). 

These nodes of WF are not displaced at change of the parameter Am, two external nodes are 
located at x-coordinates of the walls and total number of the nodes equals to m + 1. This 
property is confirmed analytically by our formula (108) and it is observed in Fig. 1 (2, 7) in [47]. 

Now let's consider, how the shape of WF of the ground state can be deformed, if in definition 
of the energy of factorization the level of the excited state at m 7^ 1 is used. In Fig. 5 the 
deformation of the shape of WF of the ground state at change of the number m of the level, 

coincident with the energy of factorization (a), and at change of the parameter A^ at m = 2 (b) 
and at m = 3 (c), is shown. In the first figure (a) the deformation of WF of the ground state 
at change of the number m is shown (such a case of deformation is absent in [47] see p. 917). 
Here, we note the following: 

• increase of the number m by 1 introduces one new perturbation (or "oscillation) into the shape 
of WF; 

• angles of leaving of WF of the ground state at points of the right and left walls keep safe at 
change of the number m, that demonstrates clearly properties 2 found in sec. 6.4. 

In the second figure (b) the deformation of WF of the ground state at change of the parameter 
A2 at m = 2 is shown. We see, that this figure after mirror refiection coincides with Fig. 1 
(6) from [47] (see p. 917), where the deformation of WF of the ground state, caused by the 
methods of the inverse problem on the basis of change of the derivative 72 of WF of the first 
excited state in the right wall, is shown. Comparing these figures, we conclude: 
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Fig. 5: The deformation of the shape of WF of the ground bound state: (a) — at change of the 
number m of the level, coincident with the energy of factorization (m = 2,3,4 at Xm = 1); (b) — at 
change of the parameter A2 (A2 = 0.01, 0.3, 1, 10, at definition of the energy of factorization the level 
with the number m = 2 is used); (c) — at change of the parameter A3 (A3 = 0.1, 0.5, 1, 10, at definition 
of the energy of factorization the level with the number m = 3 is used) 

• angles of leaving of WF of the ground state at points of the right and left walls at change of 
the parameter A2 in the method from 5.1 and at change of the coefficient 72 in the method of 
the inverse problem from [47] are not changed — this demonstrates the properties 2 found in 
sec. 6.4; 

• there is only one point, in which all deformed curves of WF of the ground state intersect between 
themselves and with the undeformed curve of this WF (without coordinates of the walls); the 
coordinate of this point is located exactly in the center of the well and coincides with the node 
of WF of the first excited state with the number m = 2, and also with zero-point of the first 
type of the deformed potential, — this visibly demonstrates properties 1 found in sec. 6.4 (this 
is the point of the first type); 

• there is the same increase of the angle between curves of the deformed WF and undeformed WF 
at point of their intersection (in the center of the well) at decreasing of the parameter A2 in the 
method from sec. 5.1 and at change of the coefficient 72 in the approach of the inverse problem 
from [47]. 

In the third figure (c) the deformation of WF of the ground state at change of the parameter 
A3 at m = 3 is shown, from here one can see, that all considered above properties of the 
deformation of WF at m = 2 are carried out and in this case (however, here WF has three 
points, in which it is not deformed: besides two points of the first type, the coordinates of 
which coincide with the nodes of WF of the excited state at the level with the number m = 3, 
which coincides with the energy of factorization, else one point of the second type appears in 
the center of the well, that confirms the properties 1 from sec. 6.4). From analysis of all figures 
one can conclude the following: 

• For arbitrary selected level with the number m, coincident with the energy of factorization, 
there are different numbers of the points of intersection between the deformed curves for the 
same WFs (using the parameter Am) for other levels with the other numbers then m (it confirms 
the properties found in sec. 6.4 and formula (112) for coordinates of these points). 
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• If the numbers m and n determine the number of oscillations of the curve of the deformed WF 
relatively its undeformed form, then using one can reinforce or enfeeble amplitude of these 
oscillations. 

The deformation of WF of the first excited state is shown in Fig. 6. In the first figure (a) 




Fig. 6: The deformation of the shape of WF of the first excited state: (a) — at change of the number 
m of the level of the deformation (m = 1,3,5, at Am = 1); (b) — at change of the parameter Am 
(A-m = 0.01,0.3,1,10, in definition of the energy of factorization the level with the number m = 1 
is used); (c) — at change of the parameter Am (Am = 0.1,0.5,1,10, in definition of the energy of 
factorization the level with the number m = 3 is used) 

the deformation of this WF at change of the number m of the level is shown (there is no such 
a case in Fig. 1 in [47], see p. 917). We see, that a behavior of this WF in such deformation 
looks like WF of the ground state (see Fig. 5 (a)): 

• increase of the number m by 1 introduce one new perturbation (or "oscillation) into the shape 
of WF; 

• angles of leaving of WF of the ground state at points of the right and left walls are not changed 
under change of the number m. 

In the second figure (b) the deformation of WF of the first excited state at change of the 
parameter Ai at m = 1 is shown. We see, that this figure at mirror reflection coincides with 
Fig. 1 (3) in [47] (see p. 917), where the deformation of WF of the first excited state, caused 
by change of the derivative 71 of WF of the ground state in the right wall using the methods 
of the inverse problem, is shown. Comparing these figures, we conclude: 

• angles of leaving of WF of the first excited state at points of the right and left walls at change 
of the parameter Ai in the method from 5.1 and at change of the coefficient 71 in the method 
of the inverse problem from [47] are not changed; 

• curves of the deformed WFs at change of Ai in the method from 5.1 and at change of 71 in the 
method of the inverse problem from [47] are not intersected (it agrees with the formula (112)). 

From the third figure (c) one can see, that the considered above properties of the deformation 
of WF of the ground state at change of Am at the selected number m, are fulfilled also for WF 
of the excited state (with increasing of m, the number of points of intersections of the deformed 
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curves of WF between themselves and with the undeformed shape of this WF for each selected 
excited state is increased, that follows also from (112)). 

In Fig. 7 one can see, which forms the deformed potential (a) and its WF of the bound state 
(b) at the level with a larger number m, coincident with the energy of factorization, have (the 
value m = 15 is chosen). 




Fig. 7: Shapes of the potential (a) and its WF of the bound state (b) at the level, coincident with 
the energy of factorization, at m = 15 and Am = 0.1 

Conclusions: 

Using the supersymmetric method from sec. 5.1, we have obtained all pictures of the de- 
formation (without displacement of the levels in the spectrum) for the rectangular well with 
infinitely high walls and its WFs of the lowest bound states, which were obtained early in the 
review [47] (see p. 916-919) on the basis of the methods of the inverse problem. Note the 
following: 

• Appropriateness of the method from sec. 5.1 for construction of new isospectral potentials 
on the basis of one given potential Vi in the form of the rectangular well with infinitely high 
external walls has proved (one can consider such a method as the method of deformation 
of the given Vi), in which for definition of the function of factorization in the construction 
of the superpotential the WF of arbitrary excited bound state (but not the ground state 
only) is used. Ah possible divergences in the found solutions can be excluded by change 
of the parameter Am inside one from regions in (117). 

• The parameter Am with arbitrary number m of the level for the bound state in the 
supersymmetric method from 5.1 plays the same role in the deformation of the potential 
and its WFs, as the derivative 7m of WF of the bound state at the same level at point 
of one from the boundaries of the potential, which is used for the deformation of the 
potential and its WFs by the methods of inverse problem (see [47] p. 917-927) (it has 
obtained at the first time). 

• Let's summarize, how WFs of the different bound states are changed under the deforma- 
tion by the parameter Am: 
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For WF at the arbitrary level with the number m, coincident with the energy of factor- 
ization: 

— WF has m + 1 nodes exactly, which are not displaced under change of and coincide 
with the nodes of the undeformed WF (96); 

— x-coordinates x = and x = a of the walls are such nodes also; 

— angles of leaving of this WF at points of the right and left walls are changed at the 
deformation. 

For WF of the other bound state at the level with arbitrary number n, not coincident 
with the energy of factorization {n ^ m): 

— There are such points, in which all deformed curves of such WF intersect between them- 
selves and with the undeformed WF (i. e. these points are not displaced at change of the 

parameter A^). This property (at the first time has been confirmed analytically by formulas 
(108) and (111)-(112), with taking into account the next point in these Conclusions) gives 
an answer on a question, stated in [47] (see p. 918): "it is not clear, why the functions 
$i(a;) and ^2{x) (i. e. WFs of the ground and the first excited states) at all values (i. e. in 
any deformation using the coefficient 72 — as analogue of A2 in the methods of the inverse 
problem) pass through only one point in the center of the well" (in the original text: " , 

$l(.x) $2(X) "). 

— One can dcvidc all these points into two types: 

* The points of the first type, x-coordinates of which equal to (113) and coincide with 
the nodes (97) or (108) of WF of the bound state at the level with the number m, 
which coincides with the energy of factorization, and also with zero-points of the first 

type of the form (101) of the potential. 

* The points of the second type, coordinates of which do not coincide with the nodes 
of WF of the bound state at the level with the number m, which coincides with 
the energy of factorization. The coordinates of these points are determined from the 
second equation of the system (112). 

— The boundary points x = and x = a are such points both the first type, and the second 
one. 

— WF of arbitrary bound state (at its deformation using the parameter A„i) has exactly m+1 
such points of the first type, and number of the points of the second type which can be 
determined from the second equation of (113). 

— Angles of leaving of such WF at points x = and x = a of the walls are not changed 
under the deformation of this WF, caused by change of the parameter A^ (or by change 

of the coefficient 7^ in the methods of the inverse problem) — in contrast to WF of the 
state with the number m\ This property (at the first time has been confirmed analytically 
in sec. 6.4) gives an answer on a question, stated in [47] (see p. 918): how to explain an 
invariable slope of WF at x-coordinates of the walls under its deformation by the coefficient 
72 using the method of the inverse problem. 

• Now let's summarize, how the shape of the well is deformed under change of the parameter 
with the fixed number m. 

According to (100), one can divide all zero-points of the potential (i. e. the points of the 
intersection of the potential with axis x) into two types (for the deformed rectangular 
well this property has been revealed and has been explained at the first time): 
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— Zero-points of the first type, which coincide with the nodes of WF of the bound state 
at the level with the number m, which coincides with the energy of factorization. These 
zero-points are not displaced Tinder the deformation of the potential by the parameter 
and their coordinates are determined from (101). 

— Zero-points of the second type, which do not coincide with the nodes of WF of the bound 
state at the level with the number m, which coincides with the energy of factorization. 
These points are displaced under the deformation of the potential by the parameter 
and their coordinates are determined from (102). 

7 Construction of new potentials, where unbound state 
is used in definition of the function of factorization 

On the example of the deformation of the rectangular well we have seen, that the supersym- 
metric method (with taking into account not only ground bound states, but excited ones also 
in construction of the superpotential) allows to obtain a whole variety of the deformations 
(without displacement of levels in spectrum), which were obtained early by the methods of the 
inverse problem. 

It turns out, that if to do not restrict oneself by the bound states only in the definition 
of the superpotential, but to construct it using a function (we shall name it as WF of the 
unbound state), which is a general solution (without implying of boundary conditions) of the 
Schrddinger equation with the starting potential Vi at arbitrary selected energy level (which can 
be not coincident with the energy levels of the potential Vi), then one can obtain new exactly 
solvable potentials (which have own energy spectra with the bound states ). Such a way for the 
construction of the superpotential enlarges further a class of exactly solvable potentials with 
only one given starting Vi (and gives more rich possibilities for the deformation of V^i), which 
was obtained early by the method from Sec. 5.1-5.2. In particular, one can construct a new 
potential, energy spectrum of which is shifted in accordence with a chosen rule relatively the 
energy spectrum of the starting potential Vi (that has not been obtained else by the other known 
SUSY QM methods). This points out to a possibility and advisability of development of SUSY 
QM methods, where in the construction of the superpotential as the energy of factorization 
the arbitrary energy (which can be not coincident with the energy spectrum of Vi) is used, 
and the function of factorization is defined concerning the unbound state at this energy. As a 
demonstration, further we shall present two new approaches of construction of new potentials, 
when as the starting potential Vi the rectangular well with infinitely high external walls in the 
form (90) is used. 

Let's consider an arbitrary unbound state at an arbitrary selected level of the energy spectrum 
of the rectangular well Vi. We denote these WF and the level by index n. One can write such 
WF in a form (at > o): 

ipi^\x) = CnSm{knX + Sn), (118) 

where kn is a wave vector, corresponding to the level {k'^ — En), Cn and 5„ are constants, 
which have complex values. 

Choosing the arbitrary energy for the definition of the energy of factorization, and WF of 
the unbound state at this energy for the definition of the function of factorization, we construct 
the superpotential (we add index w to designations of such energy, its wave vector and WF): 

W{X) = log \ip';^^\ = _l^!lliii^.i_Liy.y sin {kyjX + 6yj) = -kyj cot {kyjX + Syj) 
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or 

W{x) — —kyj cot {kyjX + Syj). (HS) 
7.1 A potential-partner to the rectangular well with a number "up" 



In the first case, we shall find a new potential V2, which is a partner to the rectangular well 
with the number "up" , obtained using direct SUSY-transformation with the superpotential in 
the form (119). 

Let's find a form of the new potential V^.: 

dW ^, d . _ 2kl 



V2{x) ^ 2—— ^ -2kyj— cot {ky,x + 5y,) - . , r \ 

dx dx sm [kyjX + dyj) 

or 

2k: 



One can see, that: 

• at Eyj > the new potential V2 has infinitely high barriers with maximums at points 

•^barrier ■ ^ ^ 

•^barrier ~ '■, ) atZ = 0,l,2... (121) 

Kyj 

• Sit Eyj < e[^^ (and wiht selected 5yj) the potential V2 becomes continuous function and 
has no any divergence in the region x e [0, a]; 

• the coordinates (121) of the divergences for V2 coincide with the coordinates of the diver- 
gences of the superpotential (119). 

However, we can take the solution (120) as the solution for the new exactly solvable potential 
V2, only if it has the bound states (without the divergences). Therefore, it is important to 
analyze WFs of V2. 

At first, let's analyze, which states for the new potential V2 ore possible at the level, coinciding 
with the energy of factorization Ey,. Direct action of operator Ay, = + W on wave function 
at the energy E^j gives zero solution. According to Sec. 4.2, there is the following non- 
normalizable solution for V2, concerned with the function ^p^^: 



^m^k±^^ /.(x)= /((^«(x'))W, (122) 



where is constant. We find a function I^: 



X X ^ \ 

Iy,(x) = J (^^^\x')fdx' ^Clj sin" {ky,x' + 5y,)dx' ^ cl j -{l - cos{2ky,x' + 25J))dx' 





Cl( , sin (2/i;^x' + 25^) 

.37 



2 \ 2/^7/ 



^ ^ 2 2kyjX + sin 25yj — sin (2kwX + 25w) 

Ik 

x'=0 ^'^'>" 
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or 



j^jjl^x; — 1^ . 

Taking into account this, we obtain: 



2/c„,a; + sin 26yj — sin {2kwX + 25^) + 



^^'^^"^^ = 4^ sin {Kx + 5^) ^^^^^ 

The solution (124) defines bound state, if the function (^^^ is continuous in the region x e [0, a] 
and equals to zero at a; = and x — a. 

Let's analyze a behavior of at point x = 0. nl we obtain: 

(2)/^\ l/l'„,Au; At, 



Akyj CI sin C^^ sin 5^ 



One can see, that at 6^, ^ wl the function 93^^ equals to zero at A^y = only. At — nl, the 
function (/?^^ (0) equals to zero also at A^ = 0. In deed, we find: 

C,„ 1 — cos 2S,„ C,i, 1 — 1 



C 2A;^( 1 — cos (2ky.x + 26^ 



We write down values 5^, and A^, when the function ip!^^ equals to zero at point x — 0: 

Xu, = 0, 6u) is arbitrary, ku, is arbitrary (/c^ 7^ 0). (125) 

Now let's analyze a behavior of the function at point x = a. If ku,a + nl, then 

from (124) at A^,, = we obtain a condition, when the function (p^^ equals to zero: 

2kyja + sin 25^ - sin (2A;^a + 25yj) = 0. (126) 

Such dependence defines a needed value of the phase 6w at selected k^, which can be arbitrary 
At k^a + Sw — ttI the function ip!^\a) equals to zero only at 

A^ = -^(2k^a + sm2S^) = 0. (127) 



In deed, we find: 



n 2k^ (1 - cos (2k^x + 25^)) 



1 - 1 Q 



2k^ {-ly 



Prom (127) we obtain the necessary condition for k^' 

sm26yj = -2ky,a. (128) 

In this case the function (f^^ equals to zero at point a; = a at some values of k-^ and Sw 

We write down such values of Sw and A^, when the function cp^^ equals to zero at points 
X = and x = a with selected k^, which can be arbitrary: 

1) /c^a + Su, ^ tt/, 2ku,a + sin 25^] - sin (2/c^a + 2Syj) = 0, A^ = 0; , . 

2) kyjU + 5yj — nl, sm25yj — —2kyja, Xyj — 0, ^ ' 

44 



where / = 0, 1, 2 . . .. 

At second, let's find WF of the arbitrary (unbound) state for V2 at the level En with the 
number n, not coinciding with the energy of factorization: 

= CnCnikn COS {knX + - ky, cot {ky,x + 5^) sin {knX + 5wi) 



or 



if^n^ {x) = CnCn (kn COS (knX + 5^) - ky, cot {ky,x + 5^) sin {knX + SrS) , (130) 

where Ay^ = + W, Cn is arbitrary constant, which we shall use instead of the factor l^^^^^ — 
E^^\~^, because we consider the non-normalizable state. 

In order to the function (130) describes the boTind state, it needs, that it must be continuous 
in the region x G [0, a] and equals to zero at points x = and x = a. Let's analyze, when this 
function can be continuous. According to (130), the divergences appear at point Xbreak at 

sin {kyjXbreak + Syj) = 0. (131) 

However, if at all points Xbreak the following condition is fulfilled: 

sin {knXbreak + ^n) = 0, (132) 

then the function (130) is continuous and has no divergences. In deed, assuming, that the 
conditions (131) and (132) are fulfilled at point x — Xbreak, we analyze the second item of the 
function (130) at this point. We obtain: 

-k cot (k X. . + 5 )sm(kx. , + 5)--k — ^^'^'^ + ^""^ '''' ^^'^'^ + ^'^^ 



sin {kwX + 5u,) 

kyj sin {kyjX + 5^,) sin {k^x + 5„) — kn cos {kyjX + 5^) cos {knX + 5^ 



kyj 

kyj cos {kyjX -\- 



7 fcn COS {kyjXbreak ~l~ ^w) COS {knXbreak ~l~ ^n) 7 /7 j_A^ ('1 V+^i- 

— I^w r 77 ~r~7 ^ — /C„ COS [KnXbreak + "nj — I, J-j ^n- 

/C^y COS \KyjXbreak ' ) 

Continuity condition of WF: the function (p!^\x), describing a state with number n and 
having the form (130), is continuous in the region x e [0,a], if each point Xbreak, defined by 
(131) at selected and S^, satisfies to the condition (132) for the corresponding kn and 5n- 
According to (130), such points Xbreak are nodes of this WF (at even /). 

Now let's find out, when the function (130) equals to zero at points x = and x = a. From 
a condition of equality to zero we obtain: 

(x = 0) kn COS Sn - kw cot 5^ sin 5„ = 0; , . 

[x = a) — >• kn cos [knO. + 5n) — kw cot {k^a + 5y,) sin (fc^a + 5^) = 0. 

Let's cos5w 7^ 0. From the first equation we obtain: 

k 

tan5„ = T^tan^tj,, cos 5„ 7^0, cos^ty^^O. (134) 
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Rewrite the second equation of the system (133) by such a way: 



kn cos 6n (cos /c„a — sin /c„a tan 5^^ — cot {k^a + 5yj) cos 5„ (sin /c„a + 7^ cos /c„a tan . 

(135) 

Prom here we obtain: 

/c„ 1 - cot (/c^a + 6^) tan 5^^, , . 

tanA;„a = -2 . (l-joj 



cot (/c^ya + ^w) + tan 



Transform this equation so: 



, 7 kn sin fc^a zio'tn 
tan/c^a = . {^^') 

Klin / , . I- s I- . K,„ 



^ COS (ky,a + Sw) cos Su, + -r^ sin (kwa + 5w) sin 5u 

k^ 



w 



In contrast to (136), numerator and denominator of the right part of the equation (137) have 
no any infinitely large value at arbitrary kn, k^ {ky, ^ 0) and Syj. Such expression is more 
convenient for analysis of zeros. 

Now let's consider a case cos (5^ = 0. From (133) wc obtain instead of the equations (134) 
an (137) the following: 

k 

tan kntt = tan k^a, cosSn — 0, cosS^j — O. (138) 

kn 

We see, that the equations (137) are (138) transcendental and in a general case it needs to 
resolve each from them numerically, that makes a problem of obtaining of new types of exactly 
solvable potentials as very difficult. In result, we shall obtain a spectrum of all possible values 
of kn (energy spectrum E^^) in dependence on selected values of k^ and 6^. However, some 
partial solutions of the new potentials (with WFs and energy spectra, in the explicit form) can 
be found more simply. Let's consider two cases. 



7.1.1 New isospectral potentials with shift of their shape along axis x 

In the first case, let's assume, that all energy levels of the new potential V2 coincide with the 
levels of the spectrum of the rectangular well Vi. Then we have: 

sinA;„a = 0. (139) 

Prom (137) at cos 5^ 7^ or from (138) at cos 5^ = we obtain: 

sinky^a^O. (140) 

We see, that the energy of factorization Ey, coincides with one level from the spectrum E^^^ for 
the rectangular well (we have an arbitrariness in a choice of the parameter 5^,). 

Now we assume, that the new potential V2, in addition to the found energy spectrum, has 
some other levels. From (138) we see, that at 6^ = ^ + ttI it is not possible (/ = 0, 1, 2 . . .). At 

7^ f + equation (137) can give such solution. As we found early, the level Eyj belongs to 
the spectrum E^\ therefore: sin/c^a — 0. Prom (137) we obtain: 

kl = -klcot^6y,. (141) 
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As for WF of the starting potential we have chosen the solutions (118) at the levels E^^ > 0, 
then k^> and A;^ > 0. Therefore, the solution (141) is possible only at complex values of 
(at ImSw 7^ 0). The complex phase gives complex values of the superpotential W and can 
lead to new types of SUSY-transformations (which one can find in literature). 

Taking into account the found energy spectrum for the potential V2, the condition of con- 
tinuity of WF of the bound state with the number n leads to the following additional inter- 
dependence between the phases S^j and 6n {w is the number of the level E^j in the spectrum 




5n^ -Syj + 27rl, 1^0,1,2... (142) 



Conclusion: At the real phase 5^, the new potential V2 has the energy spectrum E^\ coincid- 
ing completely with the energy spectrum E^^ of the rectangular well (with possible exclusion 
of one level Ey^). At Ey^ — Ei'^ (at kn — nkyj) we have the following: 

• At = ttZ (Z = 0, 1, 2 . . .) the potential V2 coincides with the potential (24) from [6] without 
shift down by value AE = Ei = |^ (see p. 278-279, Fig. 2.2; for exact comparison it needs 

to US6 Tl — 1, 2?72 — \^ L — Qj^ ^my — ^Cooper, Kharc, Sukhatmc 

+ 1, and we obtain the potential and 
spectrum without displacement: V2 = ^^cosec^(7rx/L), E^ = 2^1^ ^^)- 

— In contrast to [6] (see p. 278-279), we connect the potentials V\ and V2 by SUSY- 
transformation with non-zero lowest level £^1 7^ {Ei 7^ [6]). 

— At (5„ = ttI the condition (142) is fulfilled jointly with (134). Therefore, WF of the 
form (130) is continuous in the region x G [0, a] and equals to zero at points a; = and 
X = a. Therefore, a state described by this WF at arbitrary level with the number n (not 
coincident with Eyj) of the energy spectrum En"^ for the new potential V2 must be bound. 

— At (5„ = ttZ the conditions (1) and (2) from (129) are not fulfilled (besides = 0). 
Therefore, the state for V2 at the lowest level E^^^'^ cannot be bound and we must to 
exclude this level from the energy spectrum of the new potential V2 (this agrees with [6]). 

• At 6w ^ ttI the potential V2 displaces along axis x (using 5^) with appearance of one infinitely 
high barrier, which has a maximum with coordinate determined from (121). Here, we have two 
cases: 

— At cos 5w = the equation (142) is fulfilled jointly with the condition cos 5n = from (138) 
at 5„ = ^ + ttZ and at arbitrary odd number n. This defines a bound excited state at such 
n (continuity of WF of this state in the region x G [0, a], its zero values at points x = 

and X = a). Therefore, we obtain a new potential V2 with infinitely high barrier, maximum 
of which is located exactly in the center of the well at Xbarrier = a/2 and all levels of the 
spectrum (with the bound states) coincide with the levels with odd numbers of the spectrum 
En^ of the starting potential Vi . 

— At cos (5^ 7^ 0, for analysis of existence of excited bound states we must to analyze a 
compatibility of the conditions (142) and (134), that leads to the transcendental equation 
(134). A barrier of V2 can be displaced along axis x, its maximum has a coordinate: 



There is no bound state at the level Eyj in the potential V2. Necessary conditions (l)-(2) 
from (129) are not fulfilled. 



47 



7.1.2 New potentials with shift of the energy spectrum 

Now let's consider another case. We assume, that the levels of the new potential V2 do not 
coincide with the levels of the energy spectrum of the starting potential Vi and instead of the 
equation (139) satisfy to the following: 

coskna^O. (144) 

Let's analyze, whether the bound states in the potential V2 with such energy spectrum are 
possible. 

At cosSw 7^ from (137) with taking into account (144) we obtain: 

cos {kyjtt + 6yj) cos^u, — — -rj-sin {k^a + 6^) sin 5^. (145) 

This equation has a solution for a spectrum of values of kn with different numbers n at one 
pair of Syj and A;^, if it equals to zero. Prom here we find: 



(cos {kyjQ + Suj) = Oj ^ (sin {kwO, + S^) 7^ 0^ — >• (sin 5^ = 
(cos Syj 7^ 0^ — s> (coskyja = 0, 6yj = ttI^, 

where / = 0, 1, 2 . . .. 

At cos^u, = from (138) with taking into account (144) we obtain: 



(146) 
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COS ky^a — 0, 5yj — — + nl. (147) 



Let's write down a full solution: 



TT 

cos/c^a = 0, ^w — T^i (148) 
and the energy spectrum for the new potential V^'- 

K = ^(-1 + 2n), = kl = ^(-1 + 2n)^ (149) 

where n = 1,2, . . .. 

If at Eu, to choose the lowest level from the found spectrum then one can make sure, 
that each level of the spectrum E^^ has own bound state. A form of the potential V2 with 
such spectrum is shown in Fig. 8 (a), a shape of WF of the ground and first two excited states 
(which correspond to three lowest levels of the spectrum (149)) — in Fig. 8 (b,c). From here 
one can see, that these WF are continuous in the region x G [0, a] and equal to zero at points 
X = and x = a, that allows on their basis to define the bound states at three lowest levels for 
the potential V2. Thus, the found potential V2 in the form (120) with the parameters (I48) and 
new exactly solvable potential, which has own energy spectrum with the bound 

states. 

Note, that the spectrum of V2 is shifted relatively the spectrum of the rectangular well. 
Taking into account (93), one can write: 

_ 7rn TT _ TT 
kn — 7^ — kji, well ~ (l^'-'j 

a 2a 2a 
where kn, well is wave vector of the rectangular well with the number n. 



48 



0,2 0,4 0,6 0,8 1 

Fig. 8: The potential and its WF of the first three bound states (1=0,1,2. . . ): (a) is a shape of the 
potential at 5^ = 7r/2 + 7rl (red curve, increasing to the right) and at 5w = vr/ (blue curve, increasing to 
the left), (b) are WFs of the first three bound states (for the first three lowest levels) at 5^, = 7r/2 + 7r/, 
(c) are WFs of the first three bound states at these levels at 6u) = nl 

7.2 New potentials with insertion of one additional level in the spec- 
trum 

Now let's consider another example of the deformation of the rectangular well with infinitely 
high walls (90). We shall use the approach for the deformation of the rectangular well in Sec. 6, 
where at first we fulfill using direct SUSY-transformation a transition from the rectangular well 
(as the starting potential Vi) into a new potential V2 with transformation of one bound state 
into unbound one, and then we fulfill the inverse transition with the transformation of this 
unbound state into bound one, returning to the deformed "rectangular well" Vi(A). However, 
in contrast to the approach in Sec. 6, now in the construction of the superpotential we shall use 
an arbitrary energy Ey^ (which can be not coincident with the levels of the spectrum for Vi) 
for the definition of the energy of factorization, and the general solution of WF of the unbound 
state at this energy for the definition of the function of factorization. So, we have the function 
of factorization in the form (118) (with substitution of indexes n —>■ w) and the superpotential 
— in the form (119). 

Further, we shall restrict ourselves by finding such new potential, the energy spectrum of 
which coincides with the energy spectrum of the rectangular well with a possible insertion of 
one new level Ey^. We shall consider existence of the bound states in the found potential as a 
condition of obtaining of the new exactly solvable potential Vi(A). Therefore, first of all for the 
new potential we must to analyze a form of WF of the arbitrary state at the arbitrary level 
E!^\ which is not coincident with the energy of factorization E^. One can write such WF in 
the form (110). 

In contrast to the approach from Sec. 6, finding new solutions we shall analyze for each level 
of the starting potential Vi not WF of the bound state, but a general solution of WF of the 
unbound state (as we do it in Sec. 7), i. e. we write such WF at the level E^^ with the number 
n in the form (118). 

In result, we obtain a condition of equality of WF to zero at the level En for the new 
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potential Vi{X) at points x = and x = a (we write it without detailed computations): 



tan(5„ 



COS 2{kwa + 5w) — cos 25^, 



En ^ — E,, 



2ky^a + {sin 25^ — sin 2(/c^a + S^)^ ~^ 



:i5ii 



— Eu 



Choosing values for the parameters ky, and 5^, from (151) for each level En with the number 
n one can calculate a value of 5m providing the bound state at this level for the new potential. 
Therefore, we obtain the new potential, which has the energy spectrum, coinciding completely 
with the energy spectrum of the old potential with a possible insertion of one additional level 
Eyj, if the state at such level is bound at selected k^ and 5^,. 

For obtaining a shape of the new potential Vi{\J), the expression (99) must be modified 
with inclusion of the phase 5^: 



Vi{x, 5w) = log 2kwX + sin 25^ - sva.2{kwX + 5^,) + 2kwa\ 



(152) 



Conclusion: Inclusion of the phase 8.^, (caused by use of WF of the unbound state at the 
arbitrary energy level E^ in the construction of the superpotential) introduces an additional 
arbitrariness into the deformation of the rectangular well. 



8 Conclusions 

In paper a generalized definition of superpotential has proposed, which connects two one- 
dimensional potentials Vi and V2 with discrete energy spectra completely and where: 1) for 
definition of the energy of factorization an arbitrary level of the energy spectrum of the given 
Vi is used and the function of factorization is defined concerning a bound (ground or excited) 
state at this energy level, 2) for the definition of the energy of factorization an arbitrary energy 
(which can be not coincident with levels of the spectrum of Vi) is used and the function of 
factorization is defined concerning an unbound (or non-normalizable) state at this energy. In 
finishing, let's summarize main conclusions and perspectives. 

• It has shown, that for the unknown superpotential such its definition follows directly 
from a solution of Riccati equation at the given potential Vi, own independent variant 
of which we obtain in the paper. The found solution for the superpotential admits the 
arbitrariness on the choice of boundary conditions, imposing on wave function, which 
defines the function of factorization, and on the choice of the energy of factorization. 
This points out to appropriateness of the proposed generalization of the definition of the 
superpotential. 

• It has found, that the proposed generalization of the superpotential allows to construct 
new exactly solvable potentials with discrete energy spectra (and with own bound states) 
on the basis of one given Vi. This points out to perspectives of such generalization. 

• The generalized formalism of construction of the superpotential uses wave functions of un- 
bound and non-normalizable states. This points out perspective of further investigations 
of such states for potentials with discrete energy spectra. 
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• We improve the formalism of construction of n-parametric family of isospectral potentials, 
developed early in [6] (see p. 326-328), on the basis of the proposed generalization of the 
construction of the superpotential. 

• Using the rectangular well of finite width with infinitely high external walls as the starting 
potential Vi and constructing the superpotential, where the energy of factorization coin- 
cides with arbitrary (not only the lowest) level of the energy spectrum of the given Vi and 
the function of factorization is defined as WF of bound state at such level, by SUSY QM 
approach (the method from Sec. 5.1) we reconstruct exactly all pictures of deformation 
(without displacement of the levels) of this potential and its WFs of the lowest bound 
states, which were obtained early in review [47] (see p. 916-919) by the methods of the 
inverse problem. Note the following: 

— An appropriateness of the SUSY QM method from Sec. 5.1 for construction of 
new isospectral potentials (having own bound states without divergences) has been 
shown. 

— The parameter at arbitrary m of the level of the bound state in the SUSY QM 
method from Sec. 5.1 plays the same role in the deformation of the potential and its 

WFs, as derivative 7^ of WF of the bound state at such level in one boundary of 
the potential, used for the deformation of the potential and its WFs in the approach 
of the inverse problem (see [47] p. 917-927). 

— For WF of the bound state at the arbitrary level with the number m of the 
given Vi, if the energy of factorization £ coincides with this level {S — Em), under 
its deformation using the parameter an analysis of behavior of nodes of this WF 
is fulfilled, their coordinates are found, an analysis of angles of its leaving at points 
of right and left walls is fulfilled. 

— For WF of the bound state at the arbitrary level En with the number n, if the 
energy of factorization £ does not coincide with this level (£'„ ^ S = Em), under the 
deformation by the parameter A^ such points are found, where all deformed curves 
of this WF intersect between themselves and with the non-deformed WF, analysis 
for these points is fulfilled and their classification is proposed, it has shown, that 
angles of leaving of this WF at points a; = and not changed (in contrast 
to WF of the state with the number m, sec point above). 

— An analysis of the deformation of the shape of the well under change of the parameter 
A^ has fulfilled, an analysis of zero-points of the well (i. e. points of intersection of 
the potential with axis x) has fulfilled, their classification has proposed. 

This all makes SUSY QM methods at their efficiency practically achieved the level of the 
methods of the inverse problem. 

• Constructing the superpotential on the basis of WF of the unbound state at arbitrary 
energy of factorization (which can be located higher then the lowest level), we obtain new 
types of the deformation of the rectangular well with infinitely high external walls and 
its FWs of the lowest bound states. Using such deformation one can keep all levels of the 
energy spectrum for the given Vi at their locations (obtaining new isospectral potentials), 
or one can displace all levels of this spectrum by given rule. In particular, using only 
one superpotential, one can join two potentials, energy spectra of which (with own bound 



51 



states) have no one coincident level (it has found at the first time in the supersymmetric 
approach)! Perhaps, such a way will open a possibility to connect together boson and 
fermion components with absolutely different real energy spectra in the supersymmetric 
quantum theory of fields, using only one superpotential. 
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